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ABSTRACT 

We have observed high-dispersion echelle spectra of red giant members in the five 
open clusters NGC 1342, NGC 1662, NGC 1912, NGC 2354 and NGC 2447 and deter¬ 
mined their radial velocities and chemical compositions. These are the first chemical 
abundance measurements for all but NGC 2447. We combined our clusters from this 
and previous papers with a sample drawn from the literature for which we remeasured 
the chemical abundances to establish a common abundance scale. With this homoge¬ 
neous sample of open clusters, we study the relative elemental abundances of stars in 
open clusters in comparison with field stars as a function of age and metallicity. We 
find a range of mild enrichment of heavy (Ba—Eu) elements in young open cluster 
giants over field stars of the same metallicity. Our analysis succinct that the youngest 
stellar generations in cluster might be under-represented by the solar neighborhood 
field stars. 

Key words: - Galaxy: abundances - Galaxy: open clusters and associations - stars: 
abundances: general - open clusters: individual: NGC 1342, NGC 1662, NGC 1912, 
NGC 2354 and NGC 2447 


1 INTRODUCTION 

The first generation stars formed in galaxies are composed 
almost entirely of hydrogen and helium. When these stars 
evolve and return their nuclear-processed interiors to the 
interstellar medium (ISM), the enriched gas then present will 
be incorporated into future generation of stars. The amount 
of the chemical elements observed today and the timescales 
over which the ISM is being enriched with heavy metals are 
then, a function of many processes: the star formation rate 
(SFR), the initial mass function (IMF), the rate of element 
production and eventual return to the ISM via mass-loss 
and thresholds on the gas density for the star formation to 
proceed. The synthesis of the chemical elements and their 
return to the ISM are functions of lifetime and mass of stars. 
As the physical conditions such as the surface density of 
gas in the galactic disk and the SFR vary throughout many 
galaxies, the observed/derived abundances are a function 
of position as well. Therefore, the precise measurement of 
the variation of chemical elements as a function of radius 
in the galactic disk (i.e. the radial abundance gradient) and 
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the gradient’s temporal variation over the disk’s lifetime are 
essential to develop a complete picture of Galactic evolution. 

Stars numbering hundreds to thousands are born in 
rich open clusters (OCs) (see for example, De Silva et al. 
2009). Many high-dispersion spectroscopic abundance anal¬ 
yses support the presence of chemical homogeneity among 
cluster members with a typical star-to-star abundance scat¬ 
ter of about 0.01 to 0.05 over many elements (Carretta et 
al. 2005; De Silva et al. 2006 & 2007; Pancino et al. 2010, 
Reddy et al. 2012, hereafter Paper i and Reddy et al. 2013, 
hereafter Paper ii). This implies that the OCs harbor coeval 
group of stars all formed in a single burst of star forma¬ 
tion from a well mixed proto-cloud. An OC facilitates the 
measurement of basic parameters like age, distance, kine¬ 
matics and metallicity more accurately than for field stars. 
The broad coverage of ages and Galactocentric distances of 
these OCs make them powerful tracers to map the structure, 
kinematics, and chemistry of the Galactic disk with respect 
to Galactic coordinates and its evolution with time. 

For this reason, independent groups have attempted to 
derive two fundamental relations, the age-metallicity rela¬ 
tion and the radial metallicity gradient in the disk, using 
OC elemental abundances. But, such studies are presently 
limited by the lack of large and homogeneous datasets. Sev¬ 
eral recent attempts have been made to construct a homoge- 
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Table 1. Target clusters and their properties from the literature. 


Cluster 

1 

(deg.) 

b 

(deg.) 

Age 

(Gyr) 

[Fe/H]phot. 

(dex.) 

Rgc 

(kpc) 

(m-M)y 

(mag.) 

E(B-V) 

(mag.) 

[Fe/H],ef 

NGC 1342 

154.95 

-15.34 

0.45 

-0.16 

8.6 

10.10 

0.32 

Gratton (2000) 

NGC 1662 

187.69 

-21.11 

0.42 

-0.09 

8.4 

09.14 

0.30 

Twarog et al. (1997) 

NGC 1912 

172.25 

-1-00.69 

0.29 

-0.11 

9.1 

10.91 

0.25 

Lynga (1987) 

NGC 2354 

238.37 

-06.79 

0.13 

-0.30 

8.8 

14.01 

0.31 

Claria et al. (1999) 

NGC 2447 

240.04 

-1-00.13 

0.39 

-0.09 

8.6 

10.22 

0.05 

Claria et al. (2005) 


Table 2. The journal of the observations for each of the cluster members analysed in this paper. 


Cluster 

Star ID 

a(2000.0) 
(hh mm s) 

(5(2000.0) 

(O „ 

V 

(mag.) 

B-V 

(mag.) 

V-Ks 

(mag.) 

J-Ks 

(mag.) 

^khelio 
(km s“^) 

S/N at 
6000 A 

Date of 

observation 

NGC 1342 

4 

03 

32 

11.23 

+37 22 55.43 

09.26 

+1.31 

+3.26 

+0.78 

-10.9+0.2 

180 

29-11-2012 


6 

03 

31 

26.98 

+37 21 28.62 

09.65 

+1.18 

+3.02 

+0.66 

-10.3+0.2 

170 

29-11-2012 


7 

03 

32 

02.46 

+37 21 21.50 

09.98 

+ 1.22 

+2.89 

+0.64 

-10.8+0.2 

170 

30-11-2012 

NGC 1662 

1 

04 

48 

29.51 

+10 55 48.27 

08.34 

+1.18 

+2.99 

+0.71 

-13.6+0.2 

180 

29-11-2012 


2 

04 

48 

32.08 

+10 57 59.02 

08.87 

+1.16 

+3.02 

+0.69 

-12.9+0.2 

180 

29-11-2012 

NGC 1912 

3 

05 

28 

44.05 

+35 49 52.77 

09.85 

+1.19 

+2.85 

+0.66 

-00.2+0.2 

140 

18-11-2011 


70 

05 

29 

08.37 

+35 51 29.78 

10.04 

+1.10 

+3.13 

+0.66 

+00.6+0.2 

170 

18-11-2011 

NGC 2354 

183 

07 

13 

51.93 

-25 44 24.30 

11.41 

+1.25 

+2.79 

+0.70 

+35.6+0.4 

100 

06-03-2013 


205 

07 

13 

59.21 

-25 45 50.31 

11.13 

+1.20 

+2.50 

+0.71 

+35.0+0.3 

110 

06-03-2013 

NGC 2447 

28 

07 

44 

50.25 

-23 52 27.14 

09.96 

+0.82 

+2.33 

+0.56 

+21.1+0.3 

100 

04-03-2013 


34 

07 

44 

33.67 

-23 51 42.20 

10.15 

+0.90 

+2.21 

+0.57 

+22.7+0.3 

110 

04-03-2013 


41 

07 

44 

25.73 

-23 49 52.95 

10.16 

+0.89 

+2.28 

+0.52 

+22.2+0.2 

110 

04-03-2013 


neous sets of precise metallicity measurements for OCs (see 
for example, Heiter et al. 2014). Despite this effort to provide 
a homogeneous set of metallicities, many clusters lack spec¬ 
troscopically determined chemical compositions. This series 
is intended to examine some of these clusters. 

This is our third paper reporting a comprehensive abun¬ 
dance measurements for red giants in OCs lacking detailed 
information on their chemical composition. In our previous 
papers (Paper i & ii) we have presented abundance measure¬ 
ments for eleven OCs whose Galactocentric distances (Rgc) 
lie between 8.3 and 11.3 kpc with ages between 0.2 to 4.3 
Gyr. Here, we add five OGs, four of which have not been 
previously analysed in detail, in the Galactic anticentre di¬ 
rection. Our overall goal is to improve our understanding of 
Galactic chemical evolution (GCE). We combine our newly 
observed OGs with clusters in the literature for which high- 
resolution optical spectra have been published. Useful data 
on chemical composition is thus provided for an additional 
53 clusters. 

The layout of the paper is as follows: In Section 2 we 
describe the target selection and observations, and Section 
3 is devoted to the data reduction and radial velocity mea¬ 
surements. Section 4 is devoted to the abundance analysis 
and in Section 5 we merge our sample OCs with those from 
the literature and describe the method adopted to establish 
a common abundance scale followed by a discussion on the 
assignment of OCs to Galactic populations. We discuss in 
Section 6 the relative abundances of open clusters in com¬ 
parison with field star chemical compositions with respect 
to age and metallicity. In Section 7 we discuss the spread 
in chemical composition from cluster to cluster in the con¬ 
text of exploring the elements suitable for chemical tagging. 
Finally, Section 8 provides the conclusions. 


2 TARGET SELECTION AND OBSERVATIONS 

The new sample of OCs not yet subjected to abundance 
analysis, except NGC 2447, via high-resolution optical spec¬ 
troscopy was chosen from the New catalogue of optically vis¬ 
ible open clusters and candidate^ (Dias et al. 2002). Selec¬ 
tion of the red giant members of an OC increases the dis¬ 
tance over which the OC sample may be drawn but elimi¬ 
nates the very youngest clusters. Red giants provide spectra 
favourable for abundance determination: sharp lines with 
strengths from weak to strong for elements sampling the 
major processes of stellar nucleosynthesis. 

We have made use of the WEBDa| 3 database for the se¬ 
lection of suitable cluster members and cross-checked their 
astrometric and photometric measurements with the SIM- 
BAE0 astronomical database and more recent measurements 
are adopted. Table [1] summarizes the basic properties of the 
target clusters: references to the adopted photometric [Fe/H] 
are also given; all quantities are from the databases while 
the Galactocentric distance, Rgc, was calculated assuming 
a distance of the Sun from the Galactic centre of 8.0±0.6 
kpc (Ghez et al. 2008). The present sample of OGs covers a 
galactocentric distance of 8.4 to 9.1 kpc and an age range of 
0.1 to 0.4 Gyr. 

Observations were carried out during observing runs in 
2011 May & November, 2012 November, and 2013 March 
using the Robert G. Tull echelle spectrograph (Tull et al. 
1995) at the coude focus of the 2.7-m Harlan J. Smith tele¬ 
scope located at the McDonald observatory. On all occa- 


^ http://www.astro.iag.usp.br/ wilton/ 

^ http://www.miivle.ac.at/webda/ 

® http://simbad.u-strasbg.fr/simbad/ 
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Figure 1. Normalized spectra of red giant members of the five 
OCs as described in Table [2] are presented in descending order of 
temperature (top to bottom) in the 6120—6165 A region. All the 
spectra are on the same scale but have been displaced vertically 
from each other for clarity and better visibility. 


sions we employed a 2048x2048 24 pm pixel, backside illu¬ 
minated, anti-reflection coated CCD as a detector and the 
52.67 grooves mm“^ echelle grating with exposures centred 
at 5060 A. 

Each night’s observing rnn inclnded 5-10 zero second 
exposures (bias frames), 10-15 quartz lamp exposures (flat 
frames), and 2-3 comparison Th-Ar spectra to provide the 
reference wavelength scale. We also obtained on each ob¬ 
serving night the spectrum of hot, rapidly-rotating stars to 
monitor the presence of telluric lines. 

A total of twelve stars spread across the five OCs were 
observed. For each target, we obtained two to three expo¬ 
sures, each lasting for 20-30 min to minimize the influence 
of cosmic rays and to acquire a good signal-to-noise (S/N) 
ratio. All the spectra correspond to a resolving power of 
R > 55,000 (< 6 km s“^) as measured by the FWHM of 
Th I lines in comparison spectra, except for the members of 
NGC 2354 whose spectra was taken at the lower resolution 
of 30,000. The spectral coverage in a single exposure from 
3700 A to 9800 A across various orders was complete but 
for the inter-order gaps which begin to appear longward of 
5600 A. 


3 DATA REDUCTION AND RADIAL 
VELOCITIES 

The spectroscopic data reductions were done exactly as in 
Paper i using various routines available within the imred and 
echelle packages of the standard spectral reduction software 
iRAifl. Briefly, the two-dimensional target star frames were 

^ IRAF is a general purpose software system for the reduction 
and analysis of astronomical data distributed by NOAO, which 
is operated by the Association of Universities for Research in 


de-trended by subtracting bias level, correcting for scat¬ 
tered light and then divided by the normalized flat field. 
The individual echelle orders were traced, extracted to one¬ 
dimensional spectra and then wavelength calibrated using 
Th-Ar spectra as a reference. 

Multiple spectra were combined to acquire a single, high 
S/N spectra for each star. The combined spectra have S/N 
ratios of about 100-180 as measured around 6000 A region, 
while shortward of 5000 A the S/N ratio decreases with 
wavelength and reaches a value of about 25 around 3700 A 
region. The noisy ends of each echelle order were trimmed to 
reduce the edge effects before continuum fitting. The spec¬ 
trum was normalized interactively to unity with the avail¬ 
able cursor commands in splat task of iraf by marking con¬ 
tinuum regions, those regions which are uneffected by the 
presence of spectral lines, on each aperture which were then 
fitted by a slowly varying function such as cubic spline of 
appropriate order for better normalization. We referred our 
spectra to high-resolution spectrum of Arcturus and Sun 
(Hinkle et al. 2000) to identify the continuum regions. Spec¬ 
tra of a representative region are shown in Figure [T] for one 
star from each of the five OCs. 

We have verified the membership of target stars in each 
of the OCs using their radial velocities (RVs). RV of each 
star was calculated from its normalized spectrum by mea¬ 
suring the shift in the central wavelength from their lab¬ 
oratorary values of well defined line cores. The absence of 
systematic RV shift for the spectral lines selected from the 
blue and red wavelength regions strengthens the accuracy 
of our dispersion solution. The observed RVs were corrected 
for solar motion using iraf’s rvcorrect routine. Our mean 
RV measurements for all OCs save NGC 1912 are in fair 
agreement with the previous RV measurements for the red 
giants in OCs (Mermilliod et al. 2008). 

For NGC 1912, Mermilliod et al. (2008) derived a RV 
of —45.0±0.1 km s“^ using a single star (^405) which as 
noted by proper motion data of Mills (1967) is a cluster 
non-member. We have observed two stars (7/3 & #70) from 
this cluster, whose individual proper motions are consistent 
with the mean proper motions of cluster, implying that they 
are cluster members. Moreover, the measured RVs (Table 
[U and positions of the stars in a (B-V) colour-magnitude 
diagram (see Figure ID) are also consistent with member¬ 
ship. Based on a cross-correlation technique, Glushkova & 
Rastorguev (1991) derived a radial velocity of —1.0±0.58 
km s“^ (one star), in a good agreement with our measured 
value of (RV) = -|-0.2±0.3 km s“^ (two stars). Hence, our 
radial velocity estimate for this OC would appear to be the 
first measurement available in the literature using high dis¬ 
persion echelle spectra of two potential cluster members. 
Though the brighter giants in NGC 2447 are cluster mem¬ 
bers, we avoided them for observations as they have cool 
atmospheres with Tefj < 3800 K 

The journal of the observations for each of the cluster 
members is given in Table[2]together with the identifications, 
J2000 coordinates, available optical and 2MAS^ photom- 


Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 

® http://irsa.ipac.caltech.edu/applications/Gator 
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Figure 2. The fitting of Padova isochrones by Marigo et al. (2008) to the BV color-magnitude diagram of the clusters NGC 1342 
(photometry from Hoag et ah 1961), NGC 1662 (photometry from Hoag et al. 1961), NGC 1912 (photometry from Subramaniam & 
Sagar 1999) Sz NGC 2354 (photometry from Diirbeck 1960) using the distance, reddening and age information from Dias et al. (2002) 
catalogue. The open circles (black) represent the photometric cluster members while the filled squares (red) for the position of the 
program stars for which we derived abundances. 


etry (Cutri et al. 20030, computed heliocentric RVs, and 
S/N of the spectra around 6000 A for each of the stars. 

4 ABUNDANCE ANALYSIS 

4.1 Line selection and Equivalent widths 

The method of abundance analysis has been described in 
our Paper i, so we refer the reader to that paper for more 
information on the adopted line list and the source of refer¬ 
ences to atomic data. We employed a strict criteria to select 
clean, unblended, isolated and symmetric spectral lines due 
to weak and moderately strong lines of various atomic/ionic 
species avoiding the wavelength regions affected by telluric 
contamination and heavy line crowding. The line equiva¬ 
lent widths (EWs) were measured manually using the cursor 
commands in iraf’s splot package by fitting often a Gaus¬ 
sian profile and for a few lines a direct integration was per¬ 
formed for the best measure of EW. 

Our final line list has about 300 absorption lines cov¬ 
ering 24 elements in the spectral range 4300 — 8850 A. 


° Originally published in University of Massachusetts and In¬ 
frared Processing and Analysis Center (IPAC)/ California Insti¬ 
tute of Technology. 



Figure 3. The Htting of Padova isochrones to the BV color- 
magnitude diagram of NGC 2447 (photometry from Becker et 
al. 1976) using the distance, reddening and age information from 
Dias et al. (2002) catalogue. We marked our program stars with 
open squares (red). 
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Table 3. The linelist for all cluster members analysed in this paper. 


a(A) 

Species^ 

LEP'’ 

log gf 

#4 

NGC 1342 
#6 

#7 

NGC 

#1 

1662 

#2 

NGC 1912 
#3 #70 

NGC 2354 
#183 #205 

#28 

NGC 2447 
#34 

#41 

4668.56 

11.0 

2.10 

-1.31 

90.8 

82.3 

73.0 

79.0 

73.1 








4982.82 

11.0 

2.10 

-0.91 

116.4 

104.9 






97.1 


101.5 

101.8 

96.2 

5688.22 

11.0 

2.10 

-0.45 

149.3 

141.9 

128.8 

143.1 

134.0 

158.1 

156.5 



129.3 


130.4 

6154.23 

11.0 

2.10 

-1.55 

76.4 

66.9 

60.9 

71.9 

67.7 


96.3 

67.2 

66.0 

66.2 

62.6 

60.3 

6160.75 

11.0 

2.10 

-1.27 

96.7 

88.1 

75.4 

91.9 

83.7 

101.7 

105.8 

80.1 

87.8 

80.9 

80.5 



“ The integer part of the ’Species’ indicates the atomic number, and the decimal component indicates the ionization state 
(0 = neutral, 1 = singly ionized). 

All the lines are arranged in the order of their increasing Lower Excitation Potential (LEP). 

Note. Only a portion of this table is shown here for guidance regarding its form and content. A machine-readable version of the full 
table is available as Supporting Information with the online version of this Paper. 


Our selection criteria renders, on an average across the sam¬ 
ple of 12 stars, a list of 100 Fe i lines with lower excitation 
potentials (LEPs) ranging from 0.9 to 5.0 eV and EWs of 
up to 180 mA and 15 Fe ii lines with LEPs of 2.8 to 3.9 eV 
and EWs from ~ 25 to 130 mA. Chemical abundances for 
almost all elements are based on lines weaker than 100 mA 
and strong lines were employed only for species represented 
by a few lines (for example, Mg i). 

The accuracy of the derived abundances to a great ex¬ 
tent depend on the quality of the ^/-values which varies 
from line to line of a given element. The ^/-values of high 
accuracy (5 to 10%) are available for a large fraction of iron 
lines (Fe i & Fe ii) from recent critical reviews (Fiihr & 
Wiese 2006). This gave us a freedom to select a large sam¬ 
ple of absorption lines due to iron thanks to the presence of 
numerous iron lines in the optical region of spectra. As the 
overall goal is to have the best possible abundance deter¬ 
mination, for non-Fe elements represented by fewer lines we 
have identified a set of self-consistent absorption lines hav¬ 
ing high-quality relative <;/-values. We provide in Table O 
the atomic data adopted for each spectral line and the EW 
measurements for all stars analysed. 

We derived the solar abundances using solar EWs, mea¬ 
sured off the solar integrated disk spectrum (Kurucz et al. 
1984), and the ATLAS9 model atmosphere for Teff ,0 = 5777 
K, log go = 4.44 cm s“^ to establish a reference abundance 
scale and to reduce the systematic errors in final abun¬ 
dances. We found a microturbulence velocity of = 0.95 
km s“^ using Fe ii lines. The derived solar abundances that 
agree well with the published values of Asplund et al. (2009) 
are reported in Table 4 of Paper i. We refer to our solar abun¬ 
dances when determining the stellar abundances, [X/H] and 
[X/Fe], i.e., our analysis is essentially a differential one rel¬ 
ative to the Sun. 

4.2 Determination of atmospheric parameters 

As the EW of a spectral line is affected by the physical 
conditions in the stellar atmosphere and number density of 
absorbers, it is necessary to predetermine the atmospheric 
parameters to estimate the stellar abundances. We obtained 
our preliminary estimates from dereddenecQ optical and 

^ The adopted interstellar extinctions are (Ay, Ax, E(V-K), 
E(J-K))= (3.1, 0.28, 2.75, 0.54)’''E(B-V), where E(B-V) is taken 
from WEBDA 


2MASS photometry, following the precepts discussed in our 
Paper i. Effective temperature of each star was calculated us¬ 
ing the empirically calibrated color—temperature relation^ 
of Alonso et al. (1999). 

The surface gravities were computed by incorporating 
the known distance to the OCs, derived effective tempera¬ 
ture T^g°*, bolometric correction BCv, cluster turn-off mass 
M* and solar values of Tefi,0= 5777 K and log g0= 4.44 cm 
s“^ into the log g — TeS relation given by Allende Prieto et 
al. (1999). The BCys were derived from estimated effective 
temperatures and photometric metallicities (Table [T]) using 
Alonso et al’s (1999) calibration. 

The turn-off mass of giants has been estimated by fitting 
the cluster CMD with Padova stellar evolutionary tracks of 
Marigo et al. (2008): the adopted turn-off masses are 3.1 
Mq for NGC 1342, 2.8 M© for NGC 1662, 3.5 M© for NGC 
1912, 3.5 Mq for NGC 2354, 2.9 M© for NGC 2447. For all 
clusters, the fundamental cluster parameters (age, distance 
and reddening) of Dias et al. (2002) catalogue were adopted. 

We refined our photometric estimates of atmospheric 
parameters using spectroscopy. The line list and model at¬ 
mospheres were then used as inputs to the local thermody¬ 
namical equilibrium (LTE) line analysis and spectrum syn¬ 
thesis code MOOG0 for an abundance analysis. The one¬ 
dimensional, line-blanketed plane-parallel uniform model 
atmospheres by assumption of LTE, hydrostatic equilib¬ 
rium and flux conservation were interpolated linearly from 
the ATLAS9 model atmosphere grid of Castelli & Kurucz 
(2003). 

We performed a differential abundance analysis relative 
to the Sun by running the moog in abfind mode. Starting 
from a model with photometric estimates of atmospheric pa¬ 
rameters and the measured EWs, the individual line abun¬ 
dances were force-fitted to match the computed EWs to the 
observed ones while satisfying the following three constraints 
simultaneously: First, the microturbulence, assumed to 
be isotropic and depth independent was determined by the 
requirement that the Fe abundance from Fe ii lines be in¬ 
dependent of a line’s EW. Fe ii lines are preferred over Fe 
I as they cover a small range in LEP so that the derived 
abundance will be free of non-LTE effects, and cover a good 
range in EWs so that will be measured accurately. A check 

® Colours selected included are (B-V), (V-K) and (J-K). 

® MOOG was developed and updated by Chris Sneden and orig¬ 
inally described in Sneden (1973) 
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Table 4. Basic photometric and spectroscopic atmospheric parameters for the stars in each cluster. 


Cluster 

Star ID 

(B-V) 

:pphot 

(V-K) 

(J-K) 

1 (V-K) 

(cm s 

r|nSp©C 

(K) 

S'spec 

(cm s“^) 

^spec 

(km sec“^) 

log(L/L 0 ) 

spectroscopy 

log(L/L 0 ) 

photometry 

NGC 1342 

4 

4766 

4667 

4614 

2.21 

5100 

2.40 

1.70 

2.31 

2.34 


6 

5035 

4908 

5022 

2.50 

5200 

2.80 

1.66 

1.95 

2.14 


7 

4948 

5058 

5120 

2.68 

5200 

2.70 

1.51 

2.05 

2.03 

NGC 1662 

1 

5023 

4894 

4800 

2.32 

5100 

2.70 

1.61 

1.97 

2.28 


2 

5067 

4862 

4884 

2.53 

5200 

2.85 

1.49 

1.85 

2.06 

NGC 1912 

3 

4852 

4875 

4883 

2.27 

4950 

2.10 

1.81 

2.61 

2.42 


70 

5044 

4604 

4869 

2.28 

4950 

2.00 

1.70 

2.71 

2.30 

NGC 2354 

183 

4822 

5132 

4854 

2.43 

4920 

2.90 

1.24 

1.80 

2.35 


205 

4926 

5534 

4829 

2.47 

4850 

2.80 

1.29 

1.88 

2.44 

NGC 2447 

28 

5257 

4849 

4837 

2.58 

5050 

2.70 

1.42 

1.97 

2.02 


34 

5069 

4976 

4802 

2.64 

5050 

2.60 

1.44 

2.05 

1.98 


41 

5091 

4898 

5001 

2.64 

5100 

2.80 

1.59 

1.88 

1.97 


on the is provided by lines of Sc i, Ti i, Ti ii, V i, Cr i 
and Cr ii species. Second, the effective temperature, Teff, is 
estimated by imposing the requirement that the Fe abun¬ 
dance from Fe I lines (as they cover a good range in LEP 
0.0 to 5.0 eV) be independent of a line’s LEP (excitation 
equilibrium). This condition is also verified with the lines of 
other species like Ti i and Ni i. Third, the surface gravity 
is adjusted until the Fe abundance derived by Fe I and Fe 
II lines matched within 0.02 dex (i.e. maintaining the ion¬ 
ization equilibrium between the neutral and ionized species) 
for the derived Tefr and ^t- A check on this condition is 
also performed by Sc, Ti, V, and Cr species as they provide 
both neutral and ionized lines. As all these atmospheric pa¬ 
rameters are interdependent, several iterations are needed to 
choose a suitable model from the grid so that all the spectral 
lines in observed spectra are readily reproduced. 

The error analysis on the derived spectroscopic atmo¬ 
spheric parameters was performed in the same way as de¬ 
scribed in Paper ii. The typical errors considered in our anal¬ 
ysis are 100 K in Tefr, 0.25 cm s“^ in log g and 0.20 km s“^ 
in ^f. The derived stellar parameters for program stars in 
each of the cluster are given in Table [4] 

With few exceptions, our spectroscopic estimates are in 
good agreement with photometric ones. The uncertainties 
affecting the photometric estimates are discussed in Paper 
II, and were ascribed to adopted colors and reddening es¬ 
timates. The mean difference in photometric temperatures 
estimated using (B-V) and (V-K) is -|-53 ± 279 K and using 
(V-K) and (J-K) is -1-62 ± 236 K. The corresponding mean 
differences in and T^g”^^ with respect to 

spectroscopic Tjg^’s are -64 ± 145 K, -117 ± 290 K and 
— 179 ± 132 K respectively. Mean differences in log g and 
log(L/L©) across the 12 stars are —0.15 ± 0.22 cm-s“^ and 
+0.10 ± 0.27 cm-s“^ respectively. The corresponding com¬ 
parison of the spectroscopic [Fe/H] with photometric ones 
in Table [1] also illustrates fair agreement: A[Fe/H] = —0.02 
(NGC 1342), -0.01 (NGC 1662), 0.00 (NGG 1912), +0.11 
(NGG 2354), and -0.04 (NGG 2447). 

4.3 Abundances and error estimation 

A complete abundance analysis was conducted by running 
the abfind driver of moog with the atmospheric parameters 
determined from the Fe i and Fe ii lines as described in 


Table 5. Sensitivity of abundances to the uncertainties in the model 
parameters for the star with ID 6 in NGC 1342 with Teff= 5200 K, 
logg= 2.80 cm s“^,and ^t= 1.66 km s"^. 


Species 

Teff + 100 

K 

log g± 0. 

25 

Ct± 0.20 




^logg 



0-2 

Na I 

+0.08/ - 

0.07 

0.00/+ 0 

.03 

-0.05/+ 0.05 

0.05 

Mg I 

+0.07/ - 

0.05 

-0.01/+ 

0.05 

-0.04/+ 0.04 

0.05 

A1 I 

+0.04/ - 

0.04 

-0.01/+ 

0.01 

-0.02/+ 0.02 

0.03 

Si I 

+0.02/+ 

0.01 

+0.05/ - 

0.01 

-0.02/+ 0.02 

0.02 

Ca I 

+0.10/- 

0.09 

-0.01/+ 

0.04 

-0.02/+ 0.11 

0.07 

Sc I 

+0.11/- 

0.12 

-0.01/+ 

0.01 

-0.01/0.00 

0.07 

Sc II 

0.00/ + 0 

.01 

+0.12/ - 

0.11 

-0.06/+ 0.06 

0.07 

Ti I 

+0.12/- 

0.11 

0.00/ + 0 

.02 

-0.02/+ 0.03 

0.07 

Ti II 

+0.01/+ 

0.02 

+0.13/ - 

0.09 

-0.06/+ 0.07 

0.07 

V I 

+0.12/ - 

0.14 

-0.01/+ 

0.01 

-0.01/+ 0.01 

0.08 

Cr I 

+0.11/ - 

0.10 

+0.01/+ 

0.02 

-0.07/+ 0.08 

0.08 

Cr II 

-0.01/+ 

0.06 

+0.15/ - 

0.07 

-0.05/+ 0.07 

0.08 

Mn I 

+0.10/ - 

0.09 

+0.01/ - 

0.02 

-0.09/+ 0.11 

0.08 

Fe I 

+0.09/ - 

0.08 

+0.01/+ 

0.01 

-0.08/+ 0.10 

0.07 

Fe II 

-0.02/+ 

0.06 

+0.14/ - 

0.10 

-0.07/+ 0.08 

0.08 

Co I 

+0.08/ - 

0.08 

+0.02/ - 

0.01 

-0.03/+ 0.02 

0.05 

Ni I 

+0.06/ - 

0.06 

+0.03/ - 

0.02 

-0.05/+ 0.05 

0.05 

Cu I 

+0.08/ - 

0.06 

+0.03/ - 

0.01 

-0.05/+ 0.07 

0.05 

Zn I 

+0.01/+ 

0.02 

+0.09/ - 

0.06 

-0.13/+ 0.14 

0.09 

Y II 

+0.01/0.00 

+0.12/ - 

0.10 

-0.05/+ 0.05 

0.07 

Zr I 

+0.14/- 

0.15 

0.00/ + 0 

.01 

-0.01/+ 0.01 

0.08 

Zr II 

+0.01/+ 

0.02 

+0.13/ - 

0.09 

-0.01/+ 0.02 

0.06 

Ba II 

+0.03/ - 

0.03 

+0.08/ - 

0.08 

-0.17/ + 0.15 

0.10 

La II 

+0.03/ - 

0.01 

+0.12/ - 

0.10 

-0.02/+ 0.02 

0.07 

Ce II 

+0.02/ - 

0.01 

+0.12/ - 

0.10 

-0.05/ + 0.06 

0.07 

Nd II 

+0.03/ - 

0.01 

+0.12/ - 

0.09 

-0.01/ + 0.03 

0.06 

Sm II 

+0.03/ - 

0.01 

+0.13/ - 

0.09 

-0.02/ + 0.04 

0.07 

Eu II 

0.00/ + 0 

.01 

+0.11/ - 

0.11 

-0.02/ + 0.02 

0.06 


previous section. In most cases, the abundances are derived 
from the measured EWs but a few lines were analysed with 
synthetic spectra. 

We computed synthetic spectra for species affected by 
hyperfine structure (hfs) and isotopic splitting and/or af¬ 
fected by blends and matched them to the observed ones by 
adjustment of abundances. The features analysed by spec¬ 
trum synthesis, the adopted hfs data and isotopic ratios 
are same as those listed in our Paper ii. We have tested 
our linelists extensively to reproduce the solar and Arcturus 
spectra and measured the solar abundances to establish a 
reference abundance scale. The spectrum synthesis was car¬ 
ried out by running the MOOG in ’synth’ mode. 

Abundance results for the individual stars averaged over 
all available lines of given species are presented in Tables 
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Table 6. Elemental abundance ratios [X/Fe] for elements from Na to Eu for NGC 1342, 1662, 1912, 2354 and 2447 from this study. 
Abundances calculated by synthesis are presented in bold typeface. 


Species 

NGC 1342 

NGC 1662 

NGC 1912 

NGC 2354 

NGC 2447 

[Na I/Fe] 

-1-0.28 ±0.04 

-1-0.22 ±0.04 

-1-0.33 ±0.04 

+0.12 + 0.04 

+0.12 ± 0.04 

[Mg I/Fe] 

0.00 ± 0.04 

-0.06 ±0.03 

-1-0.03 ± 0.02 

-0.17 + 0.04 

-0.02 ± 0.04 

[A1 I/Fe] 

-0.05 ±0.02 

-0.03 ±0.03 

-1-0.06 ±0.02 

-0.11 + 0.03 

-0.14 + 0.03 

[Si I/Fe] 

-1-0.11 ±0.02 

-1-0.16 ±0.03 

-1-0.23 ±0.02 

+0.16 + 0.03 

+0.11 + 0.03 

[Ca I/Fe] 

-1-0.07 ±0.05 

-l-0.11±0.06 

-1-0.14 ±0.03 

-0.07 + 0.06 

+0.02 + 0.05 

[Sc I/Fe] 

-1-0.04 ±0.05 

-1-0.02 ±0.10 



+0.04 + 0.09 

[Sc II/Fe] 


-l-0.11±0.08 

-1-0.10 

+0.05 + 0.10 

+0.10 + 0.05 

[Ti I/Fe] 

-1-0.02 ±0.05 

-1-0.06 ±0.05 

-0.07 ±0.03 

+0.01 + 0.07 

-0.04 ± 0.05 

[Ti II/Fe] 

-0.04 ±0.05 

-1-0.05 ±0.07 

-1-0.03 ±0.06 

-0.06 + 0.06 

-0.05 + 0.06 

[V I/Fe] 

-1-0.01 ±0.06 

-1-0.03 ±0.06 

-0.07 ±0.04 

+0.04 + 0.06 

-0.02 + 0.06 

[Cr I/Fe] 

-1-0.01 ±0.06 

0.00 ± 0.04 

-1-0.01 ±0.06 

-0.03 + 0.05 

-0.04 ± 0.04 

[Cr II/Fe] 

-1-0.03 ±0.06 

-1-0.07 ±0.07 

-1-0.05 ±0.05 

-0.03 + 0.06 

+0.02 + 0.07 

[Mn I/Fe] 

-0.12 

-0.05 

-0.12 

-0.05 

-0.07 

[Fe I/H] 

-0.14 ±0.05 

-0.10 ±0.06 

-0.11 ±0.05 

-0.19 + 0.04 

-0.13 + 0.05 

[Fe II/H] 

-0.13 ±0.06 

-0.11 ±0.07 

-0.09 ±0.06 

-0.16 + 0.08 

-0.11 + 0.08 

[Co I/Fe] 

-0.03 ±0.04 

0.00 ± 0.04 

-0.10 ±0.02 

+0.07 + 0.04 

-0.04 ± 0.04 

[Ni I/Fe] 

-0.06 ±0.04 

-0.02 ±0.04 

-0.02 ±0.05 

0.00 + 0.06 

-0.07 + 0.04 

[Cu I/Fe] 

-0.29 

-0.24 

-0.30 

-0.12 

-0.28 

[Zn I/Fe] 

-0.29 

-0.13 

+0.10 

-0.31 

-0.38 

[Rb I/Fe] 

-0.04 

-0.14 

-0.30 

-0.17 

-0.18 

[Y II/Fe] 

-1-0.12 ±0.05 

-1-0.15 ±0.05 

+0.06 ±0.04 

+0.14 + 0.05 

+0.03 + 0.06 

[Zr I/Fe] 

-1-0.18 ±0.06 

-1-0.25 ±0.07 

+0.10 ±0.04 

+0.13 + 0.08 

+0.13 + 0.07 

[Zr II/Fe] 

-1-0.25 ±0.04 

-1-0.31 ±0.03 



+0.16 + 0.05 

[Ba II/Fe] 

-fO.32 

4-0.54 

+0.70 

+0.17 

+0.23 

[La II/Fe] 

-1-0.16 ±0.05 

-1-0.22 ±0.05 

+0.14 + 0.04 

+0.23 + 0.08 

+0.13 + 0.05 

[Ce II/Fe] 

-1-0.36 ±0.05 

-1-0.37 ±0.06 

+0.23 + 0.04 

+0.38 ± 0.05 

+0.32 + 0.06 

[Nd II/Fe] 

-1-0.29 ±0.04 

-1-0.26 ±0.05 

+0.13 + 0.04 

+0.33 + 0.05 

+0.22 + 0.05 

[Sm II/Fe] 

-1-0.24 ±0.05 

-1-0.22 ±0.05 

+0.04 + 0.04 

+0.24 + 0.05 

+0.19 + 0.05 

[Eu II/Fe] 

-t-0.22 

-t-0.20 

+0.07 

+0.16 

+0.22 


imi relative to solar abundances derived from the adopted 
^/-values. The Tables give the average [Fe/H] and [X/Fe] 
for all elements, standard deviation and the number of lines 
used in calculating the abundance of that element. Abun¬ 
dances calculated by synthesis are presented in bold type¬ 
face. Inspection of the Tables [9HI3 shows that, in general, 
all stars in a given cluster have very similar chemical com¬ 
positions [X/Fe] for almost all the elements and are iden¬ 
tical to within the (similar) standard deviations computed 
for an individual star. Exceptions tend to occur for species 
represented by one or a few lines, as expected when the un¬ 
certainty in measuring equivalent widths is a contributor to 
the total uncertainty. From the spread in the abundances for 
the stars of a given cluster we obtain the standard deviation 
(Ji in the Tables [SHU] in the column headed ‘average’. 


The sensitivity of the derived abundances to the vari¬ 
ations in adopted atmospheric parameters were estimated 
following the procedure described in Paper ii. The changes 
in abundances caused by varying atmospheric parameters 
by 100 K, 0.25 cm s“^ and 0.2 km s“^ with respect to the 
chosen model atmosphere are summarized in Table [S] The 
quadratic sum of all the three contributors are represented 
by CT 2 . The total error atot for each of the element is the 
quadratic sum of cti and (J 2 . The final OC mean abundances 
along with the atot from this study are presented in Table 

El 


5 COMPARISON WITH LITERATURE 

We merged our sample of sixteen OCs with the available 
high-quality results in the literature to enlarge the dataset. 
As the abundances are collected from the literature, the data 
are liable to systematic errors concerning the atomic data 
(not only the quality of 5 /-values but also the choice of spe¬ 
cific lines), model atmosphere grids employed, reference so¬ 
lar abundances adopted, abundance analysis programs and 
type of stars considered by different authors in their analy¬ 
sis. 

Many researchers have noticed a systematic offset when 
combining abundance results of OCs from different authors 
to obtain a complete picture of chemical evolution of the 
Galaxy. Among such investigators are Friel et al. (2010), 
Yong et al. (2012), and more recently the robust analysis by 
Heiter et al. (2014). In Table 13 of our Paper ii, we offered a 
comparison of our abundance estimates with literature stud¬ 
ies for the OC NGC 2682, and spotted a systematic offset 
of ±0.15 dex or smaller in [X/Fe] ratios for almost all ele¬ 
ments. Our present sample also includes a OC NGC 2447 
for which high-resolution spectroscopic abundance analyses 
have been reported by Hamdani et al. (2000) for all stars 
in common with our sample. Using their Table 4, we calcu¬ 
lated the mean abundances for all the elements and those 
values are presented in Table [7] in comparison with our re¬ 
sults. While calculating their mean abundances we have not 
included the measurement uncertainties of the order of 0.15 
dex or even larger for a few elements. We notice from the Ta- 
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Table 7. Comparison of our elemental abundance ratios [X/Fe] 
for OC NGC 2447 with those from Hamdani et al. (2000). 


Species This work Hamdani 


[Na 

i/Fe] 

-1-0.12 

± 

0.02 

+0.18 

± 

0.03 

[Mg l/Fe] 

-0.02 

± 

0.01 

+0.02 

± 

0.04 

[Si: 

[/Fe] 

-1-0.11 

± 

0.02 

0.00; 

t0.04 

[Ca 

I/Fe] 

-1-0.02 

± 

0.03 

+0.02 

± 

0.00 

[Sc 

Il/Fe] 

-1-0.10 

± 

0.01 

-0.02 

± 

0.03 

[Ti 

I/Fe] 

-0.04 

± 

0.03 

+0.09 

± 

0.03 

[V i/Fe] 

-0.02 

± 

0.02 

+0.08 

± 

0.02 

[Cr 

I/Fe] 

-0.04 

± 

0.02 

+0.07 

± 

0.01 

[Cr 

Il/Fe] 

-1-0.02 

± 

0.02 

+0.11 

± 

0.03 

[Fe 

i/H] 

-0.13 

± 

0.02 

+0.03 

± 

0.03 

[Co 

I/Fe] 

-0.04 

± 

0.03 

+0.12 

± 

0.03 

[Ni 

I/Fe] 

-0.07 

± 

0.02 

-0.11 

± 

0.02 

[Y Il/Fe] 

-1-0.03 

± 

0.02 

-0.05 

± 

0.02 

[Ce 

Il/Fe] 

-1-0.32 

± 

0.01 

+0.06 

± 

0.02 

[Eu 

Il/Fe] 

±0 

.22 

+0 

.01 


ble[7]that the mean difference in [X/Fe] between ours and 
Hamdani et al. (2000) analyses are ±0.15 dex or smaller 
for almost all elements in common: exceptions include Fe, 
Co, Ni, Ce and Eu. Therefore, we emphasize that a simple 
merger of all cluster abundances from various resources can 
easily wash out or mask subtle abundance trends unless it 
is done with extreme car 

To cancel systematic errors we have followed a different 
procedure. The only observed quantity that has to be ex¬ 
tracted from spectra for abundance estimates is the EW of 
a spectral line. On the assumption that published EWs have 
been measured reliably, we reestimated the cluster abun¬ 
dances from the EWs using our models, linelists and ref¬ 
erence solar abundances. This approach should minimize, 
or even cancel out, all sorts of systematics mentioned ear¬ 
lier and place results on a common abundance scale. A wide 
wavelength coverage of echelle data employed by various au¬ 
thors made it easy to pick up common lines for Na, Al, a- 
elements, Cr i, Fe i, Fe ii and Ni i, but note that for s- 
and r- process elements measurements are often lacking in 
the literature, thereby making it impossible to homogenize 
the data for heavy elements, but still we can measure abun¬ 
dances for [Y/Fe] for 13 OCs, [Zr/Fe] for 8 OCs and [Ce/Fe] 
for 5 OCs in the literature sample. The abundances derived 
for literature sample of OCs are shown in Tables fT51 - 1 141 

This process provides a list of 69 OCs drawn from our 
studies and from the literature, covering a range of ~ —0.5 
to 0.3 dex in metallicity and an age of few Myr to 9 Gyr. 
All these OCs, except the OC data from Villanova et al. 
(2010) and Jacobson et al. (2007 & 2011), have a high- 
quality (S/N > 50) and high-dispersion echelle spectra {R > 
25,000) observed with various telescopes and spectrographs. 
The OC data from Villanova et al. (2010), Jacobson et al. 
(2007) have high quality spectra (S/N > 50) observed at 
R ^ 17,000 and 15,000 while the OC data from Jacobson 
et al. (2011) was observed at two different resolutions R ~ 


Note, however that Mishenina et al. (2015) found a good agree¬ 
ment of their results with ours for one star in OC NGC 2506 in 
common with ours (see Table 13 in their paper), while for other 
OCs they too noticed a systematic offset in the results obtained 
by different authors. 


18,000 and 21,000. Here, we emphasize that our sample con¬ 
tribute about 23% of the total OCs explored so far with high 
quality and high-dispersion echelle spectroscopy. Moreover, 
we have done a homogeneous and comprehensive abundance 
analysis extending up to s- and r-process elements. Before 
we use this homogeneous sample to look at trends with age 
and metallicity in the disk, it is worth while to know the 
cluster’s kinematics in order to assign them to the thin or 
the thick disk. 

5.1 Assignment of OCs to Galactic populations 

The separation between thin disk, thick disk and halo pop¬ 
ulations is made by invoking the kinematic criteria used in 
studies of local field stars by Bensby et al. (2005, 2014), 
Mishenina et al. (2004), and Reddy et al. (2006). The 
method of assigning OCs either to the thin disk, the thick 
disk or the halo relies on the assumption that an observed 
sample is a mixture of three stellar populations with their 
respective Galactic space velocity UiSH, Ylsr, Y/lsr com¬ 
ponents represented by Gaussian distributions, with given 
mean values and dispersions au, <tv, <^w- The remaining 
constraints are the relative densities of thin, thick and halo 
stars in the solar vicinity. The values for the velocity dis¬ 
persions, the asymmetric drift and population fractions are 
same as those given in Table 1 of Reddy et al. (2006). 

We compute the membership probabilities Rthin, Rthick, 
Rhaio and the associated uncertainties using the equations 
and recipes described in Reddy et al.(2006). Our program 
successfully reproduces the probabilities given in Reddy et 
al. (2006). We associate the clusters to thin disk popula¬ 
tion as those that have probabilities greater than 75% and 
at least twice the probability of belonging to the thin disk 
over other populations (and likewise for other populations). 
We assume that placing a lower limit of 75% in probability 
ensures a minimum contamination of each subsample from 
the rest of the stellar populations. The clusters that do not 
satisfy these criteria but equally probable to belong to ei¬ 
ther thin and thick disks or thick disk and halo are assigned 
as respective transition objects. On account of these proba¬ 
bilities, NGC 2266 (one out of sixteen OCs studied by us) 
turned out to be a thick disk member and the remaining 
OCs in our sample display probabilities typical of the thin 
disk, not surprisingly, none belongs to the halo. Table [12] 
lists the computed space and GSR velocity components of 
OCs along with their membership probabilities. 


6 THE CLUSTER AND FIELD STAR 
POPULATIONS 

OCs form from dense molecular clouds. Stars in an OC have 
the composition of the natal cloud. It is assumed that stars 
within a cluster are of a common age and composition and 
that gas remaining with the cluster at birth is lost before it 
can be contaminated by the ejecta from the cluster’s stars. 
Such ejecta are presumed to be lost to the cluster and not 
to contaminate the remaining stars. Stars are slowly lost 
from the cluster and join the field star population. Small 
and poorly populated OCs dissolve primarily from internal 
dynamical effects in less than 10® Myr. Clusters with in¬ 
termediate masses (-^ 500 to 1000 Mq) can survive several 
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Table 8. Mean elemental abundance ratios, [X/Fe], for Na to Eu for our and literature sample of OCs, the thin disk mean abundances 
from Luck & Heiter (2007) for giants, and Luck & Heiter (2005 &; 2006) and Bensby et al. (2014) for dwarfs. Abundances calculated by 
synthesis are presented in bold typeface. 


Species 

Our sample 

Literature 
thin disk OCs 

Luck & Heiter (2007) 
thin disk giants 

Luck Heiter (2005 & 2006) 

thin disk 

Bensby et 
dwarfs 

al. (2014) 

[Pe/H] ~ 

( —U.2U to U.U) 

(-0.20 to 0.0) 

(U.U to U.2U) 

(-U.20 to U.U) 

(0.0 to 0.20) 

(-0.20 to 0.0) 

(O.U to 0.2U) 

(-U.2U to U.U) 

(U.U to U.2U) 

[Na /Fe] 

+ 0.23 ± 0.07 

+0.18+0.11 

+ 0.16+0.13 

+ 0.10 + 0.06 

+ 0.13+0.06 

+0.02+0.08 

+ 0.04+0.11 

0.00+0.05 

+ 0.02 + 0.06 

[Mg /Fe] 

+ 0.04 ± 0.08 

+0.06 + 0.10 

+ 0.03 + 0.09 

+ 0.08 + 0.10 

+ 0.08+0.10 

+0.16+0.09 

+ 0.13+0.13 

+0.07+0.11 

-0.06+0.10 

[Al /Pe] 

+ 0.03 ± 0.09 

+0.14+0.06 

+ 0.19+0.10 

+ 0.09 + 0.05 

+ 0.08+0.05 

+0.04+0.10 

+ 0.02 + 0.09 

+0.05+0.09 

+ 0.03+0.05 

[Si /Pel 

+ 0.15 ± 0.06 

+0.11+0.08 

+ 0.10+0.10 

+ 0.12 + 0.04 

+ 0.12 + 0.05 

+0.05+0.06 

+ 0.02 + 0.06 

+0.03+0.05 

+ 0.02 + 0.03 

[Ca /Fe] 

+ 0.07 ± 0.06 

+0.08+0.09 

+ 0.04+0.07 

-0.04 + 0.05 

-0.07+0.06 

+0.03+0.04 

+ 0.01 + 0.07 

+0.05+0.04 

+ 0.01 + 0.03 

[Sc /Fe] 

+ 0.07 ± 0.07 



-0.08 + 0.06 

-0.14+0.06 

+0.02+0.11 

0.00+0.09 



[Ti /Pel 

+ 0.02 ± 0.07 

+0.04+0.06 

+ 0.01 + 0.08 

0.00+0.03 

+ 0.07+0.04 

+0.05+0.07 

+ 0.05 + 0.09 

+0.04+0.06 

+ 0.01 + 0.04 

[V /Pel 

+ 0.05 ± 0.07 



-0.09 + 0.07 

-0.04+0.07 

+0.02+0.09 

0.00+0.10 



[Cr /Pe] 

+ 0.02 ± 0.05 

+0.09+0.11 

+ 0.03+0.07 

+ 0.01 + 0.05 

+ 0.03+0.03 

+0.03+0.05 

+ 0.03+0.06 

0.00+0.04 

0.00+0.02 

[Mn /Fe] 

-0.10 ± 0.07 



+ 0.06 + 0.07 

+ 0.17+0.11 

-0.05+0.12 

0.00+0.11 



[Pe /H] 

-0.12 ± 0.05 

-0.08 + 0.05 

+ 0.10+0.07 

-0.08 + 0.05 

+ 0.09+0.05 

-0.08+0.06 

+ 0.08+0.05 

-0.09+0.06 

+ 0.10+0.05 

[Co /Fe] 

+ 0.04 ± 0.07 



+ 0.06 + 0.08 

+ 0.10+0.09 

+0.03+0.10 

+ 0.04+0.07 



[Ni /Fe] 

0.00 ± 0.05 

-0.01 + 0.05 

0.00+0.06 

0.00+0.03 

+ 0.02 + 0.04 

+0.00+0.04 

+ 0.00+0.03 

-0.03+0.03 

0.00+0.04 

[Cu /Pe] 

-0.19 ± 0.07 



+ 0.01 + 0.13 

+ 0.03+0.12 

-0.12+0.19 

-0.06+0.17 



[Y /Pe] 

+ 0.09 ± 0.05 



+ 0.07 + 0.15 

+ 0.03+0.08 

-0.07+0.14 

-0.06+0.16 

-0.01 + 0.15 

-0.02+0.08 

[Ba /Fe] 

+ 0.21 ± 0.20 



+ 0.04 + 0.16 

-0.10 + 0.17 

+0.02+0.12 

-0.03+0.10 

+0.07+0.18 

+ 0.01 + 0.11 

[La /Fe] 

+ 0.17 ± 0.08 



-0.12 + 0.11 

-0.10 + 0.11 

+0.13+0.09 

+ 0.00+0.08 



[Ce /Fe] 

+ 0.22 ± 0.11 



+ 0.05 + 0.09 

-0.02 + 0.12 

+0.09+0.12 

+ 0.04+0.13 



[Nd /Fe] 

+ 0.18 ± 0.10 



-0.01 + 0.07 

-0.05 + 0.06 

+0.02+0.11 

+ 0.03+0.17 



[Sm /Fe] 

+ 0.16 ± 0.09 





-0.10+0.18 

-0.13+0.22 



[Eu /Pe] 

+ 0.13 ± O.OT 



+ 0.08 + 0.06 

+ 0.04+0.08 

+0.15+0.10 

+ 0.12 + 0.08 




Note: Although, we measured abundance ratios, [Y/Fe] (13 OCs), [Zr/Fe] (8 OCs) and [Ce/Fe] (5 OCs), for literature sample of OCs 
they were not tabulated because of the scarcity of large sample OCs with measured heavy elements. 


Gyr, if they are located in the external regions of the disk 
(Pavani & Bica 2007). One expects a close, if not exact, 
correspondence between the compositions of OCs and held 
stars. This expectation is searched for hrst through a global 
comparison of the runs of [X/Fe] with [Fe/H] and age with 
[Fe/H] for held and cluster stars. Then, in the next section, 
intracluster abundance variations, the basis of chemical tag¬ 
ging (Freeman & Bland-Hawthorn 2002), are discussed. Our 
primary sources for data on held stars are Luck & Heiter’s 
(2007) survey of local disk giants and Luck & Heiter (2005 
& 2006) and Bensby et al. (2014) for nearby FG dwarfs. 

Intrinsic variations in [X/Fe], as distinct from those aris¬ 
ing from systematic errors in abundance analyses, are in the 
main relatable to the principal sites of stellar nucleosynthesis 
which left their imprint on a cluster’s natal cloud. For the 
elements considered here, the principal sites include Type 
II supernovae from massive stars. Type la supernovae from 
white dwarfs which exceed the Chandrasekhar mass limit, 
AGB stars and possibly a merging pair of neutron stars. The 
composition of stars in a cluster will be determined by the 
cumulative effect on the natal cloud of these principal sites 
within prior generations of stars. For a cluster, the principal 
sites will seed other clouds and, thus, a future generation 
of clusters. In contrast to the case of globular clusters, it 
seems improbable that nucleosynthetically-enriched ejecta 
from cluster members will determine the composition of the 
cluster except in rare cases, e.g., the creation of Barium 
stars in a binary star with an s-process enriched AGB star 
donating mass to a less-evolved companion. 


6.1 Relative abundances and stellar ages 

A global representation of cluster and field star composi¬ 
tions is provided by Table [51 Fifteen of our 16 clusters have 
[Fe/H] from —0.2 to 0 for a mean [Fe/H] of —0.12 with the 
remaining cluster having [Fe/H] = —0.45. With the excep¬ 
tion of the latter, the clusters have been assigned to the thin 
disk. About half of the additional OGs from the literature 
have [Fe/H] ^ 0.0 and all of these belong to the thin disk. 
Table [5] provides [X/Fe] for the following nine samples: our 


15 thin disk clusters with a mean [Fe/H] = -0.1, clusters 
from the literature in the intervals —0.20 to 0.0 (15) and 0.0 
to -1-0.20 (20), thin disk giants from Luck & Heiter (2007) 
in the intervals —0.20 to 0.0 and 0.0 to -1-0.20, and thin disk 
dwarfs from both Luck & Heiter (2005 & 2006) and Bensby 
et al. (2014) for the intervals —0.20 to 0.0 and 0.0 to -1-0.20. 
Table entries come from similar but not identical analyses 
of spectra of comparable quality and, thus, small differences 
in [X/Fe] for the same element X are to be expected. Per¬ 
haps, the largest systematic differences are present between 
dwarfs and giants where the model atmospheres have ma¬ 
jor differences in Tefi and log g and, perhaps, systematic 
departures from the true atmospheres and certainly differ¬ 
ences are likely in corrections for non-LTE effects and stellar 
granulation (i.e., 3D effects). In addition, the convective en¬ 
velope of a giant brings nuclear processed material into the 
atmosphere but, except for Na, our analyses do not discuss 
elements expected to be affected such as Li, G, N, and O. 

Inspection of Table [5] shows that the [X/Fe] around the 
solar abundance [Fe/H] = 0 are similar in their mean values 
and (7 for field and OC stars for almost all elements X. For 
the purposes of discussion, elements in Table [5] are grouped 
together beginning with Na and Al and ending with the 
heavy elements Y to Eu among which are enriched in some 
OGs. The remarks below refer to thin and thick disk clusters 
and stars. 


6.2 Na and Al 

Entries for Na and Al in Table [5] are expanded in Figure (4] 
For Na, FigurejT] shows, as indicated in TableO that [Na/Fe] 
for our OCs and to a smaller degree for the literature sample 
are slightly larger than provided by Luck & Heiter (2007) 
for field giants. Both the OC and held giant samples appear 
to yield slightly more positive [Na/Fe] than the samples of 
dwarfs from Bensby et al. (2014) and Luck & Heiter (2005 
& 2006) dwarfs is not conhrmed by Bensby et al. (2014). 

That giants and dwarfs give slightly different [Na/Fe] 
values may be due to larger non-LTE effects for giants aug¬ 
mented by mild Na enrichment by the hrst dredge-up (see 
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Figure 4. The relative abundance ratios of [Na/Fe] vs. [Fe/H] 
for OCs (top panel), field giants (middle panel) and field dwarfs 
(bottom panel). Our sample of OCs are shown as green filled 
squares. Clusters from the literature are presented as red open 
squares (thin disk), blue filled circles (thick disk) and black filled 
squares (halo). Intermediate stellar populations are designated as 
magenta open circles (thin—thick disk) and cyan filled diamond 
(thick disk—halo). Luck & Heiter’s (2007) sample of field giants 
are marked as black open triangles. Bensby et al. (2014) sample of 
field dwarfs are shown as black crosses (thin disk) and red crosses 
(thick disk). Luck & Heiter (2005 & 2006) sample of field dwarfs 
are presented as blue open circles. 


also Smiljanic 2012). Non-LTE effects for the dwarfs are neg¬ 
ligible - see Bensby et al. (2014, their Figure 7). For typical 
red giants, Alexeeva et al. (2014) found non-LTE Na over¬ 
abundances from our lines to be lowered from their LTE 
values by (say) about 0.1 dex, a correction sufficient to ren¬ 
der Na in the OC giants very similar to the values seen in 
the dwarfs. Non-LTE effects for Fe i and Fe ii lines have been 
estimated by Lind et al. (2012) for giants. At the metallicity 
of our cluster’s giants (0.0 to —0.2 dex), the non-LTE correc¬ 
tions for both the Fe i and Fe ii lines are smaller than 0.02 
dex. Therefore, these corrections to our LTE Fe abundance 
do not affect the global behaviour of [Na/Fe] (or any [X/Fe]) 
with [Fe/H]. Note that rise in [Na/Fe] with [Fe/H] for [Fe/H] 
> 0.2 which occurs for both dwarfs and giants may have 
been reported first by Feltzing & Gustafsson (1998). The 
first dredge-up increases the ]Na/Fe] by a mere 0.01 dex at 
IMq, and 0.1 dex at 2Mq, to approximately 0.2 dex at 3M© 
for [Fe/H] ~ 0.0 (Karakas & Lattanzio 2014). 

In the case of Al, the samples provide a similar mean 
]A1/Fe] with the exception that the clusters drawn from the 
literature may have a higher mean [Al/Fe] (Figure [5]). This 
is not regarded as a real effect because the scatter across the 
sample is large. 

For Na and Al, our conclusion is that [Na/Fe] and 
[Al/Fe] for the thin and thick disk OCs are in fair agree¬ 
ment with respect to mean values and scatter across the 
sampled [Fe/H] range with their values for field dwarfs and 
giants. 


Figure 5. Same as Figure[4]but for [Al/Fe]. 
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Figure 6. Same as Figure!^ but for the alpha element Si. 
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Figure 7. Same as Figure [4] but for the alpha element Ca. 
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Figure 8. The relative abundance ratios [X/Fe] vs. [Fe/H] for the 
iron peak elements Cr and Ni. All the symbols have their usual 
meaning as in Figure [4] except that the OCs, field giants & field 
dwarfs are all plotted in a single panel. 


6.3 The a-elements 



mn _^_ L 
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Figure 9. The relative abundance ratios [X/Fe] versus [Fe/H] for 
the s-process element Y. All the symbols have their usual mean¬ 
ing as in Figure |4l while the red filled dots are for abundances 
of OCs from Maiorca et al. (2011) and the magenta asterisks 
are abundances of nearby young associations from D’Orazi et al. 
( 2012 ). 


Here, the a-elements comprise Mg, Si, Ca and Ti. Inspection 
of Table [H shows that the samples of cluster and field giants 
and dwarfs generally share a common [X/Fe]. A graphical 
comparison is provided in Figures|6]and[7]for two a-elements 
well represented by lines. For the giants in the field and clus¬ 
ters, [Si/Fe] appears to be about 0.1 dex greater than for the 
dwarfs. For [Ca/Fe], the cluster giants and field dwarfs share 
a very similar mean value and the sole apparent exception 
is that Luck & Heiter (2007)’s giants have a slightly lower 
value. It is of interest to note that the gradients of [Ca/Fe] 
with [Fe/H] provided by Luck & Heiter (2007)’s giants and 
Bensby et al.’s dwarfs are similar: a slight increase of [Ca/Fe] 
with decreasing [Fe/H]. [X/Fe] for a-elements in thin and 
thick disk OCs closely match their values for thin and thick 
disk field dwarfs and giants. There is a slight indication that 
the spread in the OCs [X/Fe] around [Fe/H] ~ 0.0 is larger 
than for field giants. 

A similar behaviour of a-elements Si and Ca is antici¬ 
pated because both (and Mg) are synthesized primarily in 
Type H supernovae. Although the predicted yield of Ti from 
these supernovae is too small, Ti is very often considered an 
a-element because the observed run of [Ti/Fe] with [Fe/H] 
in local thin and thick disc stars resembles those of Mg, Si 
and Ca. 


6.4 The iron group 

Within the run from Sc to Cu, Table [8] shows almost com¬ 
plete agreement for [X/Fe] from all samples. Figure [8] shows 
a typical result. Minor exceptions are noted for Mn and Cu 
where our OCs appear to have lower [X/Fe] than reported 
by Luck & Heiter for field giants and perhaps too for field 
dwarfs. A major part and even all of these exceptions may 
be attributable to neglect of hyperfine and isotopic splitting 
affecting Mn i and Cu i lines (R. E. Luck, private commu¬ 
nication). 


6.5 Heavy elements 

Even casual perusal of Table [8] shows that differences be¬ 
tween [X/Fe] for the OCs and for field stars occur for the 
heavy elements. In particular, Ba, La, Ce, Nd, Sm and pos¬ 
sibly Eu give mean [X/Fe] values that are distinctly positive 
with (TS which are greater than the corresponding as for in¬ 
dividual clusters (see, for example. Table [6]). This behaviour 
for Ba-Eu does not seem to apply to the lighter elements 
Y and Zr. This exceptional pattern for heavy elements has 
been noticed previously, see, for example, De Silva et al.’s 
(2009, Table 2) compilation of [X/Ee] and a for 24 clusters 
of which only one (M 67) is included our sample featured 
in Table |8] (Heavy element abundances were available for 
a minority of the 24 clusters.) An obvious inference is that 
OCs clearly differ in [X/Ee] for heavy elements; across our 
sample, NGC 2354 (Table fTT|l is among the most enriched in 
Ba-Sm and NGC 2682 the least enriched with the range of 
about 0.35 dex far exceeding the star-to-star a of 0.05 dex 
or less. Moreover, stars like those in NGC 2682 seem not to 
have their counterparts in the samples of field dwarfs and 
giants selected for Table [S] These field-cluster differences 
have to be related to the neutron-capture processes which 
synthesize these heavy elements. 

Heavy elements from Y to Eu sample the neutron- 
capture processes of nucleosynthesis. The weak s-process 
occuring in massive stars contributes greatly to synthesis 
of Y and Zr. The r-process in Type H SN from massive 
stars and/or in merging neutron stars is the major contrib¬ 
utor to Eu and a serious contributor to La to Sm. The main 
s-process from AGB stars is the dominant contributor to Ba 
and an appreciable contributor to La to Sm. The s-process 
contribution to solar abundances is 85% for Ba, 75% for La, 
81% for Ce, 45% for Nd and 3% for Eu (Burris et al. 2000). 

Examination of Table |8] suggests that Y has similar 
[X/Ee] in cluster and field stars (see Figures^. Although Zr 
has yet to be measured in field dwarfs and giants, it would 
appear that Zr in OCs is consistent with the claim that Y 


© 2014 RAS, MNRAS 000, [T]-?? 













12 


A. B. S. Reddy, S. Giridhar and D. L. Lambert 



Cd 

CQ 


CD 0 



Cd 

CQ 


CD 

\ 

0) 

o 


CD 


0) 

o 


0 

0 

-0 


-8 


-0 


5 
0 

5 
0 
5 
5 
0 

4 
2 
0 
2 

0 

5 

-0.4 -0.2 0 0.2 0.4 

[Fe/H] 



Figure 10. The relative abundance ratios [X/Fe] versus [Fe/H] 
for the s-process elements Ba and Ce. All the symbols have their 
usual meaning as in Figure [9] while the red filled squares are the 
Ba abundances of OCs from D’Orazi et al. (2009). 


in OCs is unchanged from its abundance in field dwarfs and 
giants, i.e., [Zr/Y] ~ 0.0. The larger scatter in [Y/Fe] from 
Luck & Heiter’s sample of dwarfs is not confirmed by Bensby 
et al. (2014). 

In contrast to Y (and Zr), the mean [X/Fe] for Ba-Sm 
for the OCs are greater than in field stars (Table [8] and 
Figure Hoi and, as we show in Section 7, there is a clear 
variation in [X/Fe] among the OCs even for the same [Fe/H]. 
Although this variation is smaller for Eu (Figure fTTI) . For 
Luck & Heiter’s sample of field dwarf stars, there may be an 
offset in [Ce/Fe] and [Eu/Fej. 


7 INTRACLUSTER ABUNDANCE 

VARIATIONS 

The proposal that Galactic field stars began life in open clus¬ 
ters implies a close match between the compositions of field 
and cluster stars. In the previous section, we showed that 
the [X/Fe] for field and cluster stars over the [Fe/H] range 
sampled by the OCs were similar with regards to mean val¬ 
ues of [X/Fe] and their a except for heavy elements from Ba 
to Eu; [X/Fe] for the latter elements span a range exceeding 
that expected purely from measurement uncertainties. For 
our OC sample, it is apparent too that the cr for [X/Fe] for 
each element is noticeably larger than the a derived from 
the several stars in a given cluster with the latter cr being 
comparable to that expected from the measurement uncer¬ 
tainties (see Table [5}. These two features suggest that there 
are measurable differences in some [X/Fe] values amongst 
clusters. 

Composition differences between clusters offer the op¬ 
portunity to exploit a technique known as ‘chemical tagging’ 
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Figure 11. Same as Figure [9] but for the r-process element Eu, 
while the red filled diamonds are for field dwarfs drawn from Koch 
& Edvardsson (2002), and the red open stars (thick disk) and 
black open stars (thin disk) are for field dwarfs from Gorbaneva 
et al. (2012). 



rCa/Fel 

Figure 12. The relative abundance ratios [Si/Fe] vs. [Ca/Fe] 
for our sample of sixteen OCs. The clusters with mean [Fe/H] 
= 0.00 ± 0.05 (one), = -0.10 ± 0.05 (eleven) and = -0.20 ± 0.05 
(three) are shown as a red filled triangle, green filled squares and 
blue filled circles respectively. The most metal-poor ([Fe/H] = 
—0.45±0.04) cluster NGC 2266 is designated as a blue open star. 

(Freeman & Bland-Hawthorn 2002). On the assumption that 
all stars in a given cluster have the same composition, the 
composition is a tag or label identifying cluster members 
even after the cluster has dissolved. The efficacy of chemi¬ 
cal tagging depends - obviously - on there being measurable 
differences in composition from cluster to cluster. 

Our sample of sixteen clusters offers an excellent op¬ 
portunity to assess the role which chemical tagging might 
play in Galactic forensics (i.e., in searching for solar twins 
and siblings, tying stars in moving groups, superclusters and 
stellar streams to potential parent open clusters). Stars in 
our sample are similar red giants analysed identically and, 
thus, systematic errors affecting the abundances and abun¬ 
dance ratios [X/Fe] should be consistent (and small) across 
the sample which spans a small range in metallicity [Fe/H], 

Composition differences which might be the basis for 
chemical tagging are most likely to be found among elements 
whose synthesis samples different nucleosynthetic processes. 
Differences should be minimized for elements sampling a sin¬ 
gle or very similar nucleosynthetic processes. Our assessment 
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[Y/Fel 

Figure 13. Same as Figure [121 but for [Si/Fe] vs. [Y/Fe]. 



rSi/Fel 

Figure 14. Same as Figure [12] but for [Nd/Fe] vs. [Si/Fe]. 

begins with elements provided by Type II and Type la super¬ 
novae, i.e., Mg to Ni. A representative illustration from this 
set of elements is in Figure [12] where [Si/Fe] versus [Ca/Fe] 
is plotted for our sample. Silicon is primarily a product of 
Type II supernovae and Fe at the typical [Fe/FI] (~ —0.1) 
of our sample is primarily a product of Type la supernovae 
and Ca is a sample of both Type II and la supernovae. For 
the subset of our OCs with [Fe/H] within —0.1 ±0.05, both 
[Si/Fe] and [Ca/Fe] cover a range of 0.13 dex when the cr for 
each cluster is ±0.03 and ±0.05 dex, respectively for Si and 
Ca. In Figure [T2| each cluster is represented by mean values 
but results for the individual stars of a cluster are closely 
arranged about the mean values for [Si/Fe] and [Ca/Fe] and 
such clusters of points are generally well separated from the 
mean values for other clusters. 

The ratio of spread to a for [Si/Fe] may suggest in¬ 
trinsic differences across the OC sample and, thus, different 
proportions of Type II to Type la products in the clusters. 
These differences are less obvious for [Ca/Fe] because of the 
larger measurement uncertainty for [Ca/Fe] and the mixed 
sites for Ca synthesis. A cautionary note is required here: 
the range in [Ni/Fe], elements with similar expected relative 
yields, is 0.17 dex for the OCs with [Fe/H] = —0.10 ± 0.05 
but a measurement a of typically ±0.03. 

Examination of Table [8] listing the [X/Fe] for field and 
cluster stars points very clearly to a difference in [X/Fe] for 
heavy elements. Among these heavy elements, the cluster- 
to-cluster scatter is small for Y and Zr, larger for Ba to 
Sm and somewhat smaller for Eu. Synthesis of Eu is pri¬ 
marily attribtable to the r-process occurring in either Type 
II supernovae and/or neutron star mergers. The r-process 
makes a contribution to the other sampled heavy elements 
(see i}631) . AGB stars and their dredge-up (and mass loss) 
are expected to be the dominant source of s-process products 


via what is known as the main s-process. Massive stars via 
the weak s-process may contribute to Y and Zr and lighter 
post-Fe elements. 

The cluster-to-cluster scatter among the Y abundances 
is not appreciably greater than that expected from measure¬ 
ment errors. In Figure fOl [Si/Fe] versus [Y/Fe] is shown for 
our clusters. With a typical error bar (see figure) of ±0.05 
and ±0.03 for Y and Si, respectively, the points in Figure 
[T3] for the 12 clusters with [Fe/H] = —0.10 ± 0.05 suggest 
a possible mild intrinsic spread in [Si/Fe] but no detectable 
intrinsic spread in [Y/Fe[. 

In contrast to the Y and Zr abundances, heavier ele¬ 
ments - Ba to Eu - show cluster-to-cluster variations larger 
than the measurement uncertainties. Figure [14] shows the 
[Nd/Fe] results as a function of [Si/Fe]. The [Nd/Fe] values 
range over 0.3 dex for our OCs with [Fe/H]= —0.10 ± 0.05 
dex with no obvious dependence on [Si/Fe]; this range is al¬ 
most ten times larger than the cr for [Nd/Fe] from a given 
cluster which is only about ±0.04 dex. 

A striking feature of the Ba to Eu variations is that 
they are not clearly related to either the s- or the r-process 
alone. Figure [TS] shows the relation between [Nd/Fe] and 
[La/Fe], [Ce/Fe[ and [Eu/Fe] for the [Fe/H]= -0.10 ± 0.05 
clusters. (Similar relations are found for [Ba/Fe] with con¬ 
siderable scatter and [Sm/Fe].) With the exception of NGC 
2682 (M 67) with an age of 4.1 Gyr, the clusters with 
[Fe/H]= —0.1 ± 0.05 have similar young ages (0.3-0.6 Gyr). 
Qualitatively, these relations suggest that either the r- and 
s-process mix for Ba-Eu across the birthplaces of these clus¬ 
ters is quite similar to a solar mix or they result from varia¬ 
tions in the Fe abundance from - presumably - Type la su¬ 
pernovae. If an s-process contribution dominated the spread, 
the [Eu/Fe] dependence on [Nd/Fe] would be much weaker; 
the s-process contribution to solar abundances is just 3% 
for Eu against 45% for Nd. Similarly, the La and Ce trends 
would be steeper with solar s-process contributions of 75% 
and 81% for La and Ce, respectively. By a similar but in¬ 
verse argument, the r-process cannot be said to dominate 
the cluster-to-cluster spread. Variations in the Fe abundance 
for a given heavy element abundance would superficially ac¬ 
count for the approximately unit slope in the three panels 
of Figure [TJ] 

Enhanced interest in heavy-element abundances in OCs 
has been stimulated by D’Orazi et al. (2009) who analysed 
by standard LTE procedures the Ba ii 5854A and 6496A 
lines in dwarf stars from 20 OCs with ages from 0.04 Gyr 
to 8.4 Gyr to find that the [Ba/Fe] values decreased with 
age: [Ba/Fe] ~ 0.0 at ages greater than about a billion years 
in OCs and field stars and rose to [Ba/Fe] ~ ±0.6 in some 
youngest clusters. 

Although the definitive study of elements across the Y- 
Eu interval for a large sample of OCs is awaited, our results 
generally confirm reports in the literature. For Y, our small 
spread in [Y/Fe] is confirmed - see, for example, Maiorca 
et al. (2011) and Mishenina et al. (2013, 2015). Figure [TS] 
shows [Y/Fe] versus age for our and other OCs. In addition, 
field giants from Luck & Heiter (2007) are included where we 
have assigned stellar ages using the relation between turn-off 
mass to obtain a red giant’s age (Takeda et al. 2008): 

log ageiyr) - 10.74-1.04 (M/M0)±O.O999 (M/M©)^ (1) 
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Figure 15. Same as Figure ^2] but for mean cluster abundances 
panel) against [Nd/Fe]. 



Age (yr) 

Figure 16. The relative abundance ratios [X/Fe] versus age (yr) 
for the s-process element Y and the r-process element Eu. All the 
symbols have their usual meaning as in Figures |9] and [TT] 

The figure shows the OC [Y/Fe] merging smoothly with 
the results for the local field giants and the spread in the two 
samples is comparable. 

Reported results for [Ba/Fe] versus [Fe/H] and age show 
a large dispersion at a given [Fe/H] or age - see, for example, 
Mishenina et al. (2013, their Figure 8) even after non-LTE 
effects on the Ba ii lines are taken into account. In contrast, 
elements such as La and Ce show tighter relations for [La/Fe] 
and [Ce/Fe] as in Figure ITTl (see also Jacobson & Friel 2013). 
It is difficult to ignore a suspicion that Ba presents a spe¬ 
cial problem because the Ba ii lines are strong whereas the 
line selection for La ii and Ce ii includes several weak lines. 
However, some authors close to the abundance analyses have 
invoked a long-forgotten neutron-capture process to account 
for the high Ba abundances. Mishenina et al. (2015) reintro¬ 
duce the intermediate neutron-capture f-process discussed 
by Cowan & Rose (1977) in which the neutron density op¬ 
erating the process is approximately 10*^® cm“®. i.e., inter¬ 
mediate between the low densities for the s-process and the 
much higher densities adopted for the r-process. However, 
abundance predictions illustrated by Mishenina at al. (2015) 
appear not to fit the OC heavy element abundances: [Ba/La] 
~ 1.2 with [Eu/La] ~ — 0.6 are predicted but with a single 
exception our results are [Ba/La] ~ 0.0 with [Eu/La] near 
- 0 . 1 . 

Intracluster abundance variations offer the opportunity 
for chemically tagging a cluster’s stars after the cluster has 


for [La/Fe] (left panel), [Ce/Fe] (middle panel) and [Eu/Fe] (right 



10 ® 10 ® 10 ^° 

Age (yr) 


Figure 17. Same as Figure but for the s-process elements La 
and Ce. 


dissolved and contributed its stars to the field population. 
As shown above, the most significant abundance variations 
are seen among the heavy elements Ba-Eu. Abundances of 
elements Na-Zn and even Y and Zr display smaller varia¬ 
tions and although such variations of some elements appear 
to exceed their measurement uncertainties, these variations 
are not obviously correlated with the tightly coupled Ba-Eu 
variations. In Figure 1181 we show ]X/Fe] for elements Na- 
Eu for the three most Ce-enriched clusters (top panel: NGC 
1662, 1342 and 2447) and the three least Ce-enriched clus¬ 
ters (bottom panel: NGC 2682, 2251 and 2482). The line 
shown in both panels connects the elements for Ce-enriched 
cluster NCC 1342. Relative to the most Ce-enriched clus¬ 
ters, the least Ce-enriched appear to be richer in Mg, Si, Cr, 
Co and Ni but by only about 0.1 dex Other elements ex¬ 
hibit even smaller differences across the spectrum of Ba-Eu 
abundance enrichment. These conclusions about intraclus¬ 
ter variations for Ba-Eu and much smaller variations for 
lighter elements confirm comments advanced by De Silva et 
al. (2009) from a cluster sample drawn from the literature. 
The present conclusions are based on our sample of OCs ob¬ 
served and analysed in a homogeneous way. Of the 16 OCs 
in the sample, the primary conclusiuons are drawn from 12 
OCs covering a narrow [Fe/H] interval (—0.10±0.05 dex), a 
limited range in Galactocentric distnce (8.3 to 9.9 kpc) and 
all with young ages but for a single cluster. 
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Figure 18. [X/Fe] versus atomic number for elements Na-Eu for the three most Ce-enriched clusters (top panel: NGC 1662, 1342 and 
2447) and the three least Ce-enriched clusters (bottom panel: NGC 2682, 2251 and 2482). The magenta line shown in both panels 
connects the elements for Ce-enriched cluster NGC 1342. 


8 CONCLUDING REMARKS 

One aim of our studies of the chemical compositions of red 
giants is to provide insights into Galactic chemical evolution 
(GCE). The principal GCE-related result presented here is 
that among our sample of young clusters with slightly sub¬ 
solar [Fe/H] there is a striking spread in the overabundances 
of elements Ba to Eu in stars with a range from [X/Fe] ~ 0.0, 
a value exhibited by held stars which are generally older, to 
approximately [X/Fe] = 0.3. (Ba is an exceptional elements 
with a higher maximum [Ba/Fe] but the great strength of 
the Ball lines may give an erroneous Ba abundance,) Such 
a spread in heavy element overabundances from cluster to 
cluster has been reported previously from young associa¬ 
tions and clusters but our analysis is the hrst to demonstrate 
clearly that the overabundances arise from contributions of 
both an s-process and an r-process. Thus, heavy elements 
from Ba to Eu are the outstanding prospects for chemical 
tagging. 

Inspection of the [X/Fe] suggests that the [X/Fe] do 
not exactly mirror the solar abundance distribution, i.e., 
the combination of s-process and r-process is not a uniform 
amplihcation of the neutron-capture processes which gave 
the solar abundances. Additionally, the fact that the high 
[X/Fe] for Ba-Eu appears in young OCs and not older OGs 
and field stars suggests that the processes are long delayed 
in making their contributions to GCE. And the existence of 
spread from OC to OC in the heavy elements [X/Fe] suggests 
that there is incomplete mixing among the OCs’ birthplaces. 
Such incompleteness may be due to the relative rarity of the 
sites for these processes and/or a long timescale for mixing. 

In the spirit of speculation to encourage further spec¬ 
troscopic analyses, the suggestion is made that the GCE is 
now influenced by r-process contributions from the mergers 
of neutron star binaries where at earlier times Type II su¬ 
pernovae were the dominant source of r-process products. 
For the s-process, if it is shown that a distinctly non-solar 
contribution is required, low mass AGB stars may now be 
more important contributors than intermediate mass AGB 
stars at earlier times. To test such speculations, it would be 
useful to determine the abundances of as many elements as 


possible for stars most enriched in Ba-Eu and then to com¬ 
pare the abundance distribution with that for stars least 
enriched in Ba-Eu. 
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Table 9. Elemental abundances for stars in the OC NGC 1342. The abundances calculated by synthesis are presented in bold typeface 
while the remaining elemental abundances were calculated using the line EWs. Numbers in the parentheses indicate the number of lines 
used in calculating the abundance of that element. 


Species 


#4 

#6 

#7 

NGC 1342Avg. 

Na I/Fe] 


+0.35 ±0.03(5) 

+0.27 + 0.03(5) 

+0.21 ± 0.03(4) 

+0.28 + 0.03 

Mg I/Fe 


±0.02 ±0.03(5) 

-0.04 + 0.05(6) 

+0.02 ± 0.05(6) 

0.00 + 0.03 

A1 I/Fe] 


-0.06 ±0.02(6) 

-0.06 + 0.01(3) 

-0.02 ± 0.03(3) 

-0.05 + 0.01 

Si I/Fe] 


±0.12 ±0.03(15) 

±0.11 ± 0.05(12) 

+0.09 + 0.04(12) 

+0.11 + 0.02 

Ca I/Fe 


±0.02 ±0.04(12) 

±0.09 + 0.03(14) 

±0.11 + 0.04(10) 

+0.07 + 0.02 

Sc I/Fe] 


-0.01 ±0.06(6) 

-0.01 + 0.07(5) 

+0.15 ± 0.05(2) 

+0.04 + 0.03 

Ti I /Fe 


-0.03 ±0.04(12) 

±0.04 + 0.04(20) 

+0.04 + 0.03(15) 

+0.02 ± 0.02 

Ti II/Fe 


-0.06 ± 0.036) 

-0.01 + 0.04(8) 

-0.06 ± 0.02(5) 

-0.04 + 0.02 

V I /Fe] 


-0.07 ±0.04(14) 

+0.04 + 0.03(13) 

+0.06 ± 0.02(9) 

+0.01 ± 0.02 

Cr I/Fe] 


-0.03 ±0.04(12) 

+0.03 ± 0.02(12) 

+0.03 ± 0.04(8) 

+0.01 ± 0.02 

Cr II/Fe 

1 

±0.02 ± 0.03(6) 

+0.04 ± 0.03(7) 

+0.04 + 0.02(6) 

+0.03 ± 0.02 

Mn I/Fe 


-0.09 

-0.10 

-0.18 

-0.12 

Fe I/H 


-0.12 ± 0.04(126) 

-0.12 + 0.03(106) 

-0.17 + 0.04(107) 

-0.14 + 0.02 

Fe II/H 


-0.10 ±0.04(19) 

-0.11 ± 0.02(19) 

-0.17 + 0.04(14) 

-0.13 ± 0.02 

Co I/Fe 


-0.06 ±0.04(9) 

0.00 + 0.05(8) 

-0.04 + 0.06(6) 

-0.03 ± 0.03 

Ni I/Fe] 


-0.10 ±0.04(33) 

-0.03 ± 0.04(32) 

-0.06 ± 0.04(32) 

-0.06 ± 0.02 

Cu I/Fe 


-0.26 

-0.30 

-0.30 

-0.29 

Zn I/Fe] 


-0.33 

-0.27 

-0.27 

-0.29 

Rb I/Fe 


-0.06 

-0.08 

±0.02 

-0.04 

Y II/Fe] 


±0.12 ±0.03(6) 

+0.13 + 0.04(8) 

+0.12 ± 0.04(6) 

+0.12 ± 0.02 

Zr I/Fe 


±0.14 ±0.04(5) 

+0.18 + 0.04(5) 

+0.21 ± 0.04(5) 

+0.18 + 0.02 

Zr II/F 

^1 

±0.24 ±0.02(3) 

+0.22 + 0.05(3) 

+0.29 ± 0.00(2) 

+0.25 ± 0.02 

Ba II/Fe 

1 

+0.36 

±0.32 

+0.27 

+0.32 

La II/Fe 


+0.17 + 0.04(4) 

+0.15 + 0.05(7) 

+0.16 ± 0.04(5) 

+0.16 ± 0.03 

Ce II/Fe 

1 

±0.30 + 0.04(5) 

+0.43 + 0.02(6) 

+0.36 ± 0.00(2) 

+0.36 ± 0.01 

Nd II/Fe] 

+0.27 + 0.03(12) 

+0.33 ± 0.02(15) 

+0.26 ± 0.03(9) 

+0.29 ± 0.02 

Sm II/F 

b1 

±0.17 + 0.03(7) 

+0.27 + 0.05(8) 

+0.27 + 0.03(7) 

+0.24 + 0.02 

Eu II/Fe 

1 

+0.21 

+0.22 

+0.22 

+0.22 


Table 10. Same as Table[9]but for stars in clusters NGC 1662 &; 1912 


Species 

NGC 1662#1 

NGC 1662#2 

NGC 1662Ayg 

1 NGC 1912#3 

NGC 1912#70 

NGC 1912Avg, 

(Na I/Fe] 

+0.23 + 0.03(4) 

+0.21 ± 0.04(4) 

+0.22 ± 0.03 

+0.33 + 0.05(4) 

+0.34 + 0.02(4) 

+0.33 + 0.04 

[Mg I/Fe] 

-0.07 + 0.02(6) 

-0.05 + 0.01(4) 

-0.06 + 0.01 

-0.01 + 0.05(4) 

+0.08 ± 0.03(5) 

+0.03 ± 0.04 

[A1 I/Fe] 

-0.02 ± 0.03(3) 

-0.05 ± 0.04(6) 

-0.03 + 0.03 

+0.05 + 0.05(4) 

+0.07 + 0.01(2) 

+0.06 + 0.04 

[Si I/Fe] 

+0.17 + 0.04(13) 

+0.15 + 0.04(12) 

+0.16 + 0.03 

+0.20 + 0.05(9) 

+0.26 + 0.05(7) 

+0.23 + 0.05 

[Ca I/Fe] 

+0.12 + 0.03(13) 

+0.11 + 0.04(13) 

+0.11 + 0.03 

+0.11 + 0.04(9) 

+0.18 + 0.04(7) 

+0.14 + 0.04 

[Sc I/Fe] 

+0.02 + 0.11(4) 


+0.02 + 0.11 




[Sc II/Fe] 

+0.15 + 0.08(5) 

+0.07 + 0.08(5) 

+0.11 + 0.06 

+0.10 

+0.10 

+0.10 

[Ti I/Fe] 

+0.07 + 0.03(16) 

+0.05 + 0.04(13) 

+0.06 + 0.03 

-0.08 + 0.05(11) 

-0.06 + 0.04(9) 

-0.07 + 0.04 

[Ti II/Fe] 

+0.03 ± 0.04(9) 

+0.06 ± 0.03(6) 

+0.05 + 0.03 

-0.07 + 0.03(3) 

+0.01 + 0.04(4) 

+0.03 + 0.03 

[V I/Fe] 

+0.04 + 0.03(11) 

+0.02 + 0.04(10) 

+0.03 + 0.03 

-0.13 + 0.03(8) 

-0.01 + 0.05(7) 

-0.07 + 0.04 

[Cr I/Fe] 

0.00 + 0.03(15) 

+0.01 + 0.03(10) 

0.00 ± 0.02 

-0.01 + 0.05(6) 

+0.04 + 0.03(4) 

+0.01 + 0.04 

[Cr II/Fe] 

+0.10 + 0.04(7) 

+0.04 + 0.04(7) 

+0.07 + 0.03 

+0.05 + 0.05(6) 

+0.06 + 0.00(2) 

+0.05 + 0.03 

[Mn I/Fe] 

-0.05 

-0.04 

-0.05 

-0.15 

-0.10 

-0.12 

[Fe I/H] 

-0.09 + 0.03(112) 

-0.11 + 0.04(107) 

-0.10 + 0.03 

-0.12 ± 0.05(49) 

-0.10 + 0.04(45) 

-0.11 + 0.04 

[Fe II/H] 

-0.10 + 0.03(15) 

-0.12 + 0.04(16) 

-0.11 ± 0.03 

-0.09 ± 0.05(13) 

-0.10 + 0.04(12) 

-0.09 + 0.04 

[Co I/Fe] 

+0.01 ± 0.04(9) 

0.00 + 0.04(8) 

0.00 + 0.03 

-0.10 + 0.01(4) 

-0.11 + 0.05(8) 

-0.10 + 0.03 

[Ni I/Fe] 

0.00 + 0.03(31) 

-0.04 + 0.04(28) 

-0.02 ± 0.03 

-0.01 ± 0.05(16) 

-0.04 + 0.04(11) 

-0.02 + 0.04 

[Cu I/Fe] 

-0.27 

-0.21 

-0.24 

-0.30 

-0.30 

-0.30 

[Zn I/Fe] 

-0.14 

-0.13 

-0.13 

+0.10 

+0.10 

+0.10 

[Rb I/Fe] 

-0.15 

-0.14 

-0.14 

-0.30 

-0.30 

-0.30 

[Y II/Fe] 

+0.18 + 0.04(6) 

+0.13 ± 0.04(5) 

+0.15 ± 0.03 

+0.04 + 0.01(4) 

+0.09 + 0.03(5) 

+0.06 + 0.02 

[Zr I/Fe] 

+0.26 ± 0.02(5) 

+0.25 ± 0.04(4) 

+0.25 ± 0.02 

+0.12 + 0.02(3) 

+0.09 + 0.05(4) 

+0.10 + 0.03 

[Zr II/Fe] 

+0.30 ± 0.04(4) 

+0.32 ± 0.05(4) 

+0.31 ± 0.03 




[Ba II/Fe] 

+0.55 

+0.54 

+0.54 

+0.70 

+0.70 

+0.70 

[La II/Fe] 

+0.23 ± 0.02(4) 

+0.21 ± 0.04(5) 

+0.22 ± 0.02 

+0.13 + 0.02(4) 

+0.16 + 0.03(3) 

+0.14 + 0.02 

[Ce II/Fe] 

+0.37 + 0.04(5) 

+0.38 ± 0.02(6) 

+0.37 + 0.02 

+0.24 + 0.02(3) 

+0.23 + 0.03(3) 

+0.23 + 0.02 

[Nd II/Fe] 

+0.25 ± 0.03(13) 

+0.28 + 0.03(13) 

+0.26 ± 0.02 

+0.14 + 0.03(6) 

+0.13 + 0.04(9) 

+0.13 + 0.03 

[Sm II/Fe] 

+0.24 + 0.03(7) 

+0.20 + 0.04(8) 

+0.22 ± 0.03 

+0.06 + 0.04(4) 

+0.03 + 0.03(5) 

+0.04 + 0.03 

[Eu II/Fe] 

+0.20 

+0.20 

+0.20 

+0.10 

+0.05 

+0.07 
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Table 11. Same as Table[9]but for stars in the OC NGC 2354 


Species 

NGC 2354#183 

NGC 2354#205 

NGC 2354Avg. 


Na I/Fe] 

+0.14 ±0.04(3) 

±0.11 + 0.03(2) 

±0.12 + 0.03 


Mg I/Fe] 

-0.18 ±0.04(3) 

-0.16 + 0.04(5) 

-0.17 + 0.03 


A1 I/Fe] 

-0.14 ±0.05(5) 

-0.09 + 0.03(3) 

-0.11 + 0.03 


Si I/Fe] 

±0.15 ±0.03(8) 

±0.17 + 0.03(7) 

±0.16 + 0.02 


Ca I/Fe] 

-0.07 ±0.05(9) 

-0.07 + 0.04(5) 

-0.07 + 0.03 


Sc II/Fe] 

±0.05 ±0.10(3) 


±0.05 + 0.10 


Ti I/Fe] 

-0.07 ±0.04(12) 

±0.09 + 0.04(16) 

±0.01 + 0.03 


Ti II/Fe] 

-0.12 ±0.03(3) 

0.00 ± 0.03(4) 

-0.06 + 0.02 


V I/Fe] 

0.00 ±0.05(8) 

±0.08 + 0.05(7) 

±0.04 ± 0.04 


Cr I/Fe] 

-0.07 ±0.04(5) 

±0.01 + 0.04(6) 

-0.03 + 0.03 


Cr II/Fe] 

-0.03 ±0.03(2) 

-0.03 + 0.04(2) 

-0.03 + 0.03 


Mn I/Fe] 

-0.01 

-0.09 

-0.05 


Fe I/H] 

-0.20 ±0.03(71) 

-0.18 + 0.03(67) 

-0.19 + 0.02 


Fe II/H] 

-0.18 ±0.02(5) 

-0.15 + 0.03(8) 

-0.16 + 0.02 


Co I/Fe] 

0.00 ±0.03(5) 

±0.14 + 0.01(3) 

±0.07 + 0.02 


Ni I/Fe] 

-0.01 ±0.04(18) 

±0.01 + 0.04(13) 

0.00 ± 0.03 


Cu I/Fe] 

-0.13 

-0.11 

-0.12 


Zn I/Fe] 

-0.30 

-0.33 

-0.31 


Rb I/Fe] 

-0.21 

-0.14 

-0.17 


Y II/Fe] 

±0.16 ±0.03(3) 

±0.12 + 0.03(3) 

±0.14 + 0.02 


Zr I/Fe] 

±0.11 ±0.04(3) 

±0.16 + 0.01(3) 

±0.13 + 0.02 


Ba II/Fe] 

+0.19 

±0.16 

+0.17 


La II/Fe] 

+0.25 ± 0.04(4) 

±0.22 ± 0.04(4) 

±0.23 + 0.03 


Ce II/Fe] 

+0.34 + 0.03(4) 

±0.43 ± 0.03(4) 

±0.38 + 0.02 


Nd II/Fe] 

±0.35 + 0.04(4) 

±0.32 + 0.02(8) 

±0.33 + 0.02 


Sm II/Fe] 

±0.26 + 0.02(3) 

±0.22 + 0.02(6) 

±0.24 + 0.01 


Eu II/Fe] 


±0.16 

+0.16 


Table 12. Same as Table [9] but for stars in the OC NGC 2447. 


Species 


NGC 2447#28 

NGC 2447#34 

NGC 2447#41 

NGC 2447Avg. 

]Na I/Fe] 


±0.13 + 0.04(5) 

+0.13 + 0.04(3) 

+0.09 ± 0.02(4) 

+0.12 + 0.02 

Mg I/Fe 


-0.08 + 0.02(5) 

+0.01 + 0.02(4) 

0.00 ± 0.02(4) 

-0.02 + 0.01 

A1 I/Fe] 


-0.16 ±0.03^) 

-0.13 +0.03b) 

-0.12 + 0.03(6) 

-0.14 + 0.02 

Si I/Fe] 


+0.09 ± 0.05(8) 

+0.11 ± 0.04b) 

+0.12 + 0.03(12) 

+0.11 + 0.02 

Ca I/Fe] 


±0.02 + 0.04(10) 

+0.01 ± 0.04b) 

+0.02 + 0.05^0) 

+0.02 ± 0.03 

Sc I/Fe] 




+0.04 + 0.09(2) 

+0.04 ± 0.09 

Sc II/Fe 




+o.io±o.oib) 

+0.10 + 0.01 

Ti I/Fe] 


-0.05 + 0.04(19) 

-0.06 + 0.04(17) 

-0.01 + 0.05(17) 

-0.04 + 0.03 

Ti II/Fe 


-0.05 ± 0.045) 

-0.05 + 0.04(8) 

-0.06 + 0.03(6) 

-0.05 + 0.02 

V I/Fe] 


-0.03 + 0.04(13) 

-0.04 + 0.04(14) 

+0.02 + 0.03(16) 

-0.02 + 0.02 

Cr I/Fe] 


-0.05 + 0.04(13) 

-0.01 + 0.03^0) 

-0.06 + 0.03(9) 

-0.04 ± 0.02 

Cr II/Fe 


+0.03 ± 0.02(3) 

0.00 + 0.02(5) 

+0.02 ± 0.04b) 

+0.02 ± 0.02 

Mn I/Fe 


-0.05 

-0.08 

-0.07 

-0.07 

Fe I/H] 


-0.11 + 0.04(101) 

-0.13 + 0.04(102) 

-0.14 + 0.04(120) 

-0.13 + 0.02 

Fe II/H] 


-0.09 ± 0.04(9) 

-0.12 + 0.03(11) 

-0.13 + 0.04(16) 

-0.11 + 0.02 

Co I/Fe] 


-0.07 + 0.04(6) 

-0.01 + 0.05(7) 

-0.03 ± 0.04(6) 

-0.04 ± 0.03 

Ni I/Fe] 


-0.08 + 0.04(25) 

-0.07 + 0.04(30) 

-0.05 ± 0.04(27) 

-0.07 + 0.02 

Cu I/Fe] 


-0.27 

-0.25 

-0.32 

-0.28 

Zn I/Fe] 


-0.34 

-0.32 

-0.49 

-0.38 

Rb I/Fe] 


-0.20 

-0.18 

-0.17 

-0.18 

Y II/Fe] 


-0.01 + 0.04(7) 

-0.02 + 0.03(5) 

+0.12 + 0.05(7) 

+0.03 ± 0.02 

Zr I/Fe] 


+0.12 +0.02b) 

+0.11 ± 0.05b) 

+0.16 +0.04b) 

+0.13 + 0.02 

Zr II/Fe 



+0.12 +0.04b) 

+0.20±0.00b) 

+0.16 + 0.03 

Ba II/Fe 




+0.23 

+0.23 

La II/Fe 


+0.15 + 0.04(5) 

+0.11 + 0.05(3) 

+0.12 + 0.03(6) 

+0.13 + 0.02 

Ce II/Fe 


+0.32±0.0lb) 

+0.32 +0.03b) 

+0.31 +0.03b) 

+0.32 + 0.01 

Nd II/Fe] 

+0.22 + 0.05(14) 

+0.22 + 0.05(10) 

+0.22 + 0.04(14) 

+0.22 + 0.03 

Sm II/Fe 

1 

+0.22 ± 0.05(6) 

+0.16 + 0.05(5) 

+0.20 ± 0.04(5) 

+0.19 + 0.03 

[Eu II/Fe] 

+0.20 

+0.22 

+0.23 

+0.22 


© 2014 RAS, MNRAS 000, [T]-?? 
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Table 13. Elemental abundance ratios [X/Fe] for elements Na, Al, Mg, Si, Ca, Ti, Cr and Ni for the literature sample. 


Cluster 

[Na/Fe] 

[Mg/Fe] 

[Al/Fe] 

[Si/Fe] 

[Ca/Fe] 

[Ti/Fe] 

[Cr/Fe] 

[Fe/H] 

[Ni/Fe] 

NGC 6404 


0.08±0.07 

0.27±0.07 

Thin disk 

0.18±0.10 -0.11±0.10 

-0.07±0.07 

0.05±0.07 

0.03±0.10 

0.02±0.07 

NGC 6583 


0.03±0.01 

0.22±0.01 

0.01±0.13 

O.OlibO.OO 

-0.04±0.02 

-0.09±0.01 

0.30±0.01 

0.03±0.04 

NGC 3960 

0.06±0.02 

-0.06±0.04 

-0.02±0.02 

0.02±0.02 

0.03±0.03 

0.01±0.03 

0.02±0.02 

0.12±0.02 

-0.05±0.02 

Collinder 261 

0.08±0.06 

0.17±0.07 

0.33±0.06 

0.24±0.05 

-0.01±0.06 

-0.02±0.09 

-0.01±0.08 

0.07±0.06 

0.05±0.06 

NGC 6192 


0.04±0.03 

0.13±0.03 

0.11±0.07 

0.01±0.05 

0.03±0.09 

0.00±0.08 

0.09±0.09 

-0.05±0.07 

IC 4756 

0.18±0.08 

0.00±0.03 


0.00±0.06 

0.04±0.06 

-0.07±0.05 

0.03±0.04 

0.04±0.04 

-0.07±0.04 

NGC 3532 

0.27±0.06 

0.02±0.01 


0.02±0.05 

0.08±0.05 

-0.07±0.05 

-0.01±0.04 

0.05±0.04 

-0.07±0.03 

NGC 6281 

0.23±0.01 

0.05±0.02 


0.07±0.02 

0.10±0.05 

-0.07±0.04 

-0.02±0.06 

0.05±0.05 

-0.02±0.07 

NGC 6633 

0.24±0.01 

0.05±0.00 


0.06±0.03 

O.OlibO.06 

-0.07±0.08 

0.05±0.07 

0.00±0.07 

-0.03±0.05 

IC 2602 

0.06±0.01 



0.06±0.04 

0.08±0.06 

0.04±0.00 

0.06±0.00 

-0.06±0.05 

-0.01ib0.06 

IC 2391 

0.03±0.03 



0.04±0.02 

0.07±0.03 

0.10±0.06 

0.02±0.00 

-0.02±0.03 

0.03±0.02 

NGC 7160 

-0.20±0.08 



0.07±0.09 

-0.02±0.09 



0.13±0.06 

0.01±0.06 

NGC 5822 

0.23±0.04 


-0.Olio.00 

0.07±0.01 

0.03±0.08 

-0.02±0.13 

0.02io.ll 

0.01±0.09 

-0.04±0.07 

NGC 6134 

0.26±0.08 

-0.03±0.05 


0.02±0.02 

0.04±0.09 

0.04±0.02 

0.10±0.03 

0.15±0.04 

0.01±0.06 

NGC 6253 

0.29±0.25 

-0.05±0.12 

0.30ib.l6 

0.23±0.08 

-0.02±0.02 

0.07±0.08 

0.11±0.07 

0.18±0.01 

0.08±0.03 

NGC 6791 

0.15±0.10 



0.23±0.05 

0.11±0.12 

0.26±0.08 

0.10±0.06 

0.24±0.10 

0.18±0.03 

NGC 6819 

0.43±0.01 

-0.Olio.01 

0.15i0.09 

0.18±0.09 

0.19±0.08 

0.07±0.03 

0.18±0.03 

0.04±0.09 

-0.04±0.01 

NGC 2660 

0.07±0.03 

0.00±0.03 

O.lOibO.06 

0.03±0.01 

0.02±0.04 

0.07±0.02 

0.02±0.02 

0.06±0.02 

-0.04±0.02 

NGC 188 

0.20±0.06 

0.25±0.04 

0.25±0.06 

0.24±0.04 

0.02±0.02 

0.07±0.06 


0.11±0.07 

0.04±0.09 

NGC 2477 

0.15±0.02 

-0.04±0.03 

0.12±0.03 

0.11±0.04 

0.05±0.03 

0.01±0.06 

0.05i0.06 

0.07±0.04 

-0.03±0.05 

NGC 6939 

0.12±0.02 


0.27±0.01 

0.10±0.04 

0.16±0.03 



0.07±0.01 

0.02±0.02 

Collinder 110 

0.13±0.03 

0.04io.01 

0.20±0.01 

0.05±0.03 

0.03±0.02 

0.08±0.03 

0.08i0.03 

-0.07±0.02 

-0.03±0.04 

NGC 2099 

0.19±0.00 

-0.10±0.08 

0.13±0.07 

0.02±0.00 

0.00±0.02 

-0.02±0.04 

0.00±0.08 

-0.03±0.02 

-0.01±0.02 

NGC 2420 

0.09±0.06 

-0.11±0.04 

0.11±0.03 

-0.02±0.04 

0.02±0.01 

0.12±0.08 

0.00±0.08 

-0.06±0.02 

O.OlihO.Ol 

IC 2714 

0.23±0.03 

0.04±0.02 


0.22±0.08 

0.07±0.05 

0.01±0.03 

-0.06±0.07 

0.08±0.07 

-0.07±0.04 

NGC 7789 

0.07±0.05 

-0.09±0.03 

0.12io.06 

-0.12±0.03 

-0.05±0.03 

-0.04±0.04 

-0.01±0.03 

0.06±0.09 

0.03±0.09 

NGC 1245 

0.18±0.05 

0.16±0.06 


0.10±0.04 

0.09±0.04 



0.03±0.05 

0.00±0.04 

NGC 2194 

0.16±0.03 

0.11±0.08 


0.12±0.02 

O.OlibO.06 



-0.10±0.03 

-0.03±0.06 

NGC 2355 

0.21±0.06 

0.27±0.03 

o.iiio.oo 

0.19±0.04 

0.23±0.10 



-0.12±0.01 

0.03±0.09 

NGC 7142 

0.50±0.09 


0.27±0.12 

0.26±0.04 

0.20±0.06 



-0.03ih0.ll 

0.07±0.08 

NGC 2324 

0.25±0.04 

0.15io.04 

0.04±0.05 

0.08±0.05 

0.11±0.07 

0.02±0.09 

0.05i0.03 

-0.12±0.04 

-0.05±0.03 

NGC 2141 

0.26±0.12 

0.08±0.10 

0.09±0.01 

0.20±0.16 

0.09±0.15 

0.07±0.10 


-0.14±0.16 

-0.08ih0.ll 

NGC 2158 

0.20±0.08 

0.09±0.11 

0.16±0.10 

0.20±0.22 

0.29±0.08 

0.05±0.08 


-0.15±0.18 

0.05±0.21 

Berkeley 75 

0.28±0.04 




0.16±0.06 

O.llitO.ll 

-o.ioio.io 

-0.28±0.18 

0.04±0.06 

Berkeley 39 

0.14±0.06 

0.12io.06 

0.2lio.03 

0.20±0.04 

0.04±0.07 

0.13±0.02 

0.36±0.05 

-0.15±0.09 

0.02±0.09 

NGC 1883 

O.lOitO.ll 

0.05±0.05 

0.06±0.14 

0.16±0.02 

-0.01±0.03 

0.01±0.02 

0.24±0.07 

-0.06±0.01 

-0.07±0.04 

Ruprecht 4 

O.lOitO.Ol 

-0.04±0.00 

0.13±0.00 

-0.01±0.08 

O.llitO.ll 

0.02±0.06 

0.00±0.13 

-0.04±0.05 

-0.09±0.08 

Melotte 66 

0.10±0.09 

0.12±0.03 

0.32±0.04 

Thin - Thick disk 

0.20±0.02 0.07±0.02 

0.02±0.02 

-0.01±0.02 

-0.16±0.01 

0.00±0.02 

NGC 2425 

0.09±0.07 

-0.03±0.07 


0.02±0.07 

0.18±0.03 



O.lOihO.Ol 

-0.05±0.06 

Tombaugh 2 




Thick disk 

0.20±0.11 



-0.31±0.10 

0.06±0.12 

Berkeley 73 

0.19±0.06 

O.OOitO.OO 

0.17±0.14 

0.14±0.07 

-0.07±0.07 

0.03±0.07 

0.05±0.13 

-0.23±0.10 

0.04±0.14 

Berkeley 32 

-0.06±0.04 

0.13±0.05 

0.11±0.06 

0.10±0.04 

0.07±0.06 

0.07±0.05 

0.01±0.04 

-0.22±0.05 

0.01±0.04 

NGC 1193 

0.14±0.02 

0.25±0.12 

0.13±0.14 

0.17±0.07 

0.03±0.12 



-0.18±0.10 

0.07±0.11 

Berkeley 22 

0.14±0.07 

0.12±0.05 

0.41±0.06 

0.16±0.09 

-0.02±0.03 

0.17±0.07 


-0.35±0.02 

0.01±0.04 

Berkeley 18 

0.12±0.03 

0.20±0.06 

0.30±0.02 

0.11±0.02 

0.07±0.02 

0.12±0.00 


-0.35±0.00 

-0.02±0.05 

NGC 2243 

0.12±0.05 

0.22±0.04 

0.40±0.09 

0.19±0.11 

0.13±0.03 



-0.42±0.09 

0.02±0.09 

Berkeley 20 

0.12±0.05 

0.26±0.03 

0.28±0.02 

0.12±0.10 

0.05±0.07 

-0.01±0.09 

-0.03i0.05 

-0.27±0.05 

0.03±0.11 

Berkeley 21 

0.18±0.06 

0.10±0.08 

0.34±0.02 

0.18±0.03 

0.05±0.10 

0.05±0.09 


-0.32±0.04 

0.03±0.06 

Saurer 1 

0.21±0.07 


0.41±0.04 

0.26±0.01 

0.25±0.04 

0.25±0.05 


-0.39±0.01 

0.20±0.01 

Berkeley 33 

0.17±0.02 

-0.12±0.10 

0.33±0.00 

Thick disk 

0.10±0.06 

- Halo 

0.12±0.13 

-0.10±0.09 

0.01±0.12 

-0.34±0.09 

-0.12±0.08 

Berkeley 25 

0.08±0.14 

O.lOitO.OO 


Halo 

0.22±0.17 

0.05±0.10 

0.20±0.13 

0.22±0.14 

-0.30±0.05 

-0.01±0.21 

Berkeley 29 

0.11±0.05 

0.18±0.03 

0.23i0.10 

0.11±0.06 

0.06±0.04 

0.10±0.12 

-0.09±0.05 

-0.37±0.05 

0.00±0.04 

Berkeley 31 

0.14±0.02 

0.11±0.03 

0.18±0.05 

0.16±0.09 

0.08±0.03 

-0.01±0.05 


-0.24±0.06 

0.01±0.02 


Table 14. Elemental abundance ratios [X/Fe] for elements Y, Zr and Ce for the literature sample. 


Cluster 

[Y/Fe] 

[Zr/Fe] 

[Ce/Fe] 

Cluster 

[Y/Fe] 

[Zr/Fe] [Ce/Fe] 

IC 4756 

Thin disk 

0.06±0.09 

0.23±0.06 


Thick disk 

NGC 3532 

0.04±0.11 


0.29±0.08 

Berkeley 32 


-0.21±0.07 

NGC 6281 

0.05±0.08 


0.26±0.02 

Berkeley 22 


-0.44±0.01 

NGC 6633 

0.15±0.11 


0.23ih0.09 

Berkeley 18 


-0.28ih0.02 

NGC 5822 

0.02±0.11 


0.20±0.08 

Berkeley 21 


-0.26±0.10 

NGC 6134 

0.23±0.01 



Thick disk 

- Halo 

NGC 6791 

0.02±0.03 



Berkeley 33 

-0.22ih0.05 


NGC 188 
Collinder 110 
NGC 2099 
NGC 2420 
NGC 7789 
NGC 2141 

NGC 2158 
NGC 1883 
Ruprecht 4 

0.00i0.02 

0.02±0.05 

-0.04ih0.08 

0.27±0.02 

-0.05i0.08 

0.08io.ll 

-0.05i0.01 

-0.12±0.04 

-0.08±0.05 
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Table 15. The input data used in calculating the OC’s space coordinates and velocity components. The full sample contains 69 OCs 
with the parameters extracted from DAML catalogue, while those from different sources in the literature are shown in bold numbers. 


Cluster 

a(2000.0) 

hh:mm:ss 

5(2000.0) 

O ' " 

1 

(deg.) 

b 

(deg.) 

d© 

(kpc) 

RV 

(km sec”^) 

/iQ COS 5 
(mas yr~^) 

(mas yr~ ^) 

Rgc 

(kpc) 

Age 

(Myr) 

NGC 752 

01:57:41 

+37:47:06 

137.125 

-23.254 

0.457+0.091 

6.30+0.10 

7.50+0.32 

-11.50+0.32 

8.3 

1122 

NGC 1817 

05:12:15 

+ 16:41:24 

186.156 

-13.096 

1.972+0.394 

66.40+0.11 

1.02+1.06 

-6.51 + 1.06 

9.9 

409 

NGC 2360 

07:17:43 

-15:38:30 

229.807 

-1.426 

1.887+0.377 

29.55+0.30 

-3.51 + 1.49 

8.07+1.57 

9.3 

561 

NGC 2506 

08:00:01 

-10:46:12 

230.564 

9.935 

3.460+0.692 

84.63+0.18 

-2.55+0.20 

0.37+0.14 

10.5 

1109 

NGC 2266 

06:43:19 

+26:58:12 

187.790 

10.294 

3.400+0.680 

-29.70+0.20 

-1.98+0.59 

-4.32+0.59 

11.3 

1200 

NGC 2335 

07:06:49 

-10:01:42 

223.600 

-1.183 

1.417+0.283 

-3.21+0.10 

-0.91 + 1.56 

-3.18+1.25 

9.1 

162 

NGC 2482 

07:55:12 

-24:15:30 

241.626 

2.035 

1.343+0.269 

39.00+0.20 

-4.93+0.53 

1.60+0.53 

8.7 

402 

NGC 2251 

06:34:38 

+08:22:00 

203.58 

+00.10 

1.349+0.270 

26.25+0.14 

-0.24+0.55 

-2.60+0.51 

9.2 

267 

NGC 2527 

08:04:58 

-28:08:48 

246.087 

1.855 

0.601+0.120 

40.55+0.14 

-4.10+0.55 

6.40+0.55 

8.3 

446 

NGC 2539 

08:10:37 

-12:49:06 

233.705 

11.111 

1.363+0.273 

29.20+0.14 

-4.07+0.27 

-1.83+0.27 

8.9 

371 

NGC 2682 

08:51:18 

+ 11:48:00 

215.696 

31.896 

0.908+0.182 

34.83+0.10 

-8.62+0.28 

-6.00+0.28 

8.7 

4300 

NGC 1342 

03:31:38 

+37:22:36 

154.952 

-15.342 

0.665+0.133 

-10.67+0.11 

-1.15+0.87 

-2.80+0.87 

8.6 

452 

NGC 1662 

04:48:27 

+ 10:56:12 

187.695 

-21.114 

0.437+0.087 

-13.35+0.29 

-1.93+0.28 

-2.10+0.28 

8.4 

422 

NGC 2447 

07:44:30 

-23:51:24 

240.039 

0.135 

1.037+0.207 

22.08+0.18 

-4.85+0.33 

4.47+0.33 

8.6 

387 

NGC 2354 

07:14:10 

-25:41:24 

238.368 

-6.792 

1.445+0.289 

34.40+0.26 

-2.44+0.64 

-2.96+0.64 

8.8 

134 

NGC 1912 

05:28:40 

+35:50:54 

172.250 

0.695 

1.066+0.213 

0.18+0.19 

0.23+0.17 

-5.44+0.19 

9.1 

290 

NGC 1193 

03:05:56 

+44:23:00 

146.812 

-12.163 

4.571+0.914 

-88.10+0.20 

1.41+0.54 

-4.04+0.54 

12.0 

5012 

NGC 1245 

03:14:42 

+47:14:12 

146.647 

-8.931 

3 . 000 + 0 . 300 ^ 

- 29 . 70 + 1 . 20 ^ 

-2.98+0.41 

-3.05+0.36 

10.6 

1047 

NGC 188 

04:72:08 

+85:15:18 

122.843 

22.384 

1.700+0.340 

-42.36+0.04 

-1.48+1.25 

-0.56+1.24 

9.0 

4285 

Berkeley 18 

05:22:12 

+45:24:00 

163.632 

5.017 

5.200+1.040 

-5.50+1.10 

-2.63+0.35 

-5.95+0.35 

13.1 

4266 

NGC 1883 

05:25:54 

+46:29:24 

163.083 

6.159 

4.800+0.960 

-30.80+0.60 

-1.62+0.39 

-1.23+0.39 

12.7 

1000 

Berkeley 20 

05:32:37 

+00:11:18 

203.483 

-17.373 

8.400+1.680 

70.00+13.00 

1.51+0.81 

-4.11+0.81 

15.9 

6026 

Berkeley 21 

05:51:42 

+21:47:00 

186.840 

-2.509 

6.200+1.240 

-1.00+1.00 

-0.65+0.53 

-6.20+0.53 

14.2 

2188 

NGC 2099 

05:52:18 

+32:33:12 

177.635 

3.091 

1.383+0.277 

8.30+0.20 

3.78+0.29 

-7.09+0.29 

9.4 

347 

Berkeley 22 

05:58:27 

+07:45:24 

199.877 

-8.079 

6 .000+1.200 

95.30+0.20 

- 0 . 40 + 3 . 87 ^ 

- 3 . 81 + 2 . 90 “* 

13.8 

3311 

NGC 2141 

06:02:55 

+ 10:26:48 

198.044 

-5.811 

4.033+0.807 

24.10+0.30 

1.38+0.38 

-1.78+0.38 

11.9 

1702 

NGC 2158 

06:07:25 

+24:05:48 

186.634 

1.781 

4 . 036 + 0 . 125 ^ 

26.90+1.90 

1.43+0.29 

-3.28+0.20 

12.0 

1054 

NGC 2194 

06:13:45 

+ 12:48:24 

197.250 

-2.349 

1 . 903 + 0 . 120 ^ 

7 . 80 + 0 . 80 ^ 

-0.31+0.64 

-4.39+0.64 

9.8 

327 

Berkeley 73 

06:22:00 

-06:21:00 

215.278 

-9.424 

9.800+1.960 

95.70+0.20 

+ 0 . 54 + 2 . 67 '* 

- 2 . 60 + 2 . 91 “* 

16.9 

1514 

NGC 2243 

06:29:34 

-31:17:00 

239.478 

-18.014 

4.458+0.892 

59.84+0.41 

2.53+0.54 

2.90+1.30 

10.9 

1076 

Collinder 110 

06:38:24 

+02:01:00 

120.996 

8.602 

2.184+0.437 

41 . 00 + 3 . 80 ^ 

0.59+0.39 

-1.41+0.39 

9.3 

1412 

Berkeley 25 

06:41:00 

-16:31:00 

226.612 

-9.700 

11.400+2.280 

134.30+0.20 

- 1 . 27 + 4 . 13 ^ 

- 2 . 53 + 3 . 03 “* 

17.8 

5012 

Ruprecht 4 

06:48:54 

-10:32:00 

222.047 

-5.339 

4.700+0.940 

47.50+1.00 

- 0 . 33 + 2 . 34 ^ 

+ 0 . 93 + 2 . 64 ^ 

11.9 

794 

Berkeley 75 

06:49:01 

-23:59:48 

234.307 

-11.188 

9.100+1.820 

94.60+0.35 

- 0 . 61 + 3 . 37 “* 

+ 0 . 56 + 3 . 36 “* 

15.2 

3981 

Berkeley 29 

06:53:18 

+ 16:55:00 

197.983 

8.025 

14.871+2.974 

24.80+0.10 

-0.14+0.80 

-4.75+0.58 

22.6 

1059 

Berkeley 31 

06:57:36 

+08:16:00 

206.254 

5.120 

8.272+1.654 

55.70+0.70 

-4.30+0.52 

-3.97+0.52 

15.8 

2056 

Berkeley 32 

06:58:06 

+06:26:00 

207.952 

4.404 

3 . 420 + 0 . 684 ^ 

105.20+0.86 

-2.78+0.88 

-3.21+0.88 

11.1 

3388 

Tombaugh 2 

07:03:05 

-20:49:00 

232.832 

-6.880 

6.080+1.216 

120.90+0.40 

- 2 . 89 + 2 . 65 “* 

0 . 95 + 2 . 07 ^ 

12.6 

1995 

NGC 2324 

07:04:07 

+01:02:42 

213.447 

3.297 

3.800+0.760 

41.81+0.22 

-1.70+1.01 

-2.77+1.01 

11.4 

631 

NGC 2355 

07:16:59 

+ 13:45:00 

203.390 

11.803 

1 . 928 + 0 . 130 ^ 

35.02+0.16 

-2.50+0.80 

-3.00+1.50 

9.8 

162 

Saurer 1 

07:20:56 

+01:48:29 

214.689 

7.386 

13.200+2.640 

98.00+9.00 

+ 0 . 48 + 2 . 89 “* 

- 0 . 69 + 3 . 55 “* 

20.3 

5012 

Melotte 66 

07:26:23 

-47:40:00 

259.559 

-14.244 

4.313+0.863 

23.00+6.00 

-4.18+0.61 

7.67+1.56 

9.7 

2786 

NGC 2425 

07:38:17 

-14:52:42 

231.504 

3.297 

3 . 062 + 0 . 100 ^ 

102.90±1.20 

-3.64+0.87 

-0.76+0.87 

10.2 

2188 

NGC 2420 

07:38:23 

+21:34:24 

198.107 

19.634 

2.480+0.496 

73.57+0.15 

-1.32+0.42 

-4.18+0.26 

10.3 

1995 

Berkeley 33 

06:57:42 

-13:13:00 

225.424 

-4.622 

6500+1.300 

76.60+0.50 

-5.73+0.96 

+3.71+0.96 

13.4 

800 

Berkeley 39 

07:46:42 

-04:36:00 

223.462 

10.095 

4.780+0.956 

55.00+7.00 

- 1 . 62 + 3 . 61 “* 

- 1 . 66 + 3 . 07 “* 

11.9 

7943 

NGC 2477 

07:52:10 

-38:31:48 

253.563 

-5.838 

1.300+0.260 

7.26+0.12 

-0.73+0.57 

1.22+0.54 

8.5 

603 

NGC 2660 

08:42:38 

-47:12:00 

265.929 

-3.013 

2.826+0.565 

21.34+0.46 

-5.82+0.81 

7.40+0.83 

8.7 

1079 

IC 2602 

10:42:58 

-64:24:00 

289.601 

-4.906 

0.161+0.032 

10 .10+0.22 

-17.49+0.22 

10 .10+0.22 

7.9 

32 

NGC 3532 

11:05:39 

-58:45:12 

289.571 

1.347 

0.492+0.098 

4.33+0.34 

-10.95+0.28 

4.80+0.28 

7.8 

300 

IC 2714 

11:17:27 

-62:44:00 

198.107 

19.634 

1.320+0.264 

-14.10+0.30 

-2.56+0.73 

-3.40+0.61 

9.2 

348 

NGC 3960 

11:50:33 

-55:40:24 

294.367 

6.183 

1.850+0.370 

-22.26+0.36 

-7.01+0.24 

-0.45+0.33 

7.4 

1259 

IC 2391 

08:40:32 

-53:02:00 

270.362 

-06.839 

0.175+0.035 

14.40+0.10 

-33.21+0.30 

-5.97+0.30 

8.0 

46 

Collinder 261 

12:37:57 

-68:22:00 

301.684 

-5.528 

2.190+0.438 

-30.00+9.00 

-6.53+0.72 

-1.04+0.72 

7.1 

8912 

NGC 5822 

15:04:21 

-54:23:48 

321.573 

3.593 

0.933+0.187 

-29.31+0.18 

-7.95+0.24 

-8.20+0.24 

7.3 

891 

NGC 6134 

16:27:46 

-49:09:06 

334.917 

-0.198 

0.913+0.183 

-25.70+0.19 

-0.86+0.88 

-4.60+0.88 

7.2 

929 

NGC 6192 

16:40:23 

-43:22:00 

340.647 

2.122 

1.547+0.309 

-7.93+0.21 

3.73+0.83 

3.18+1.39 

6.6 

135 

NGC 6253 

16:59:05 

-52:42:30 

335.460 

-6.251 

1.510+0.302 

-29.40+1.30 

- 1 . 55 + 2 . 27 “* 

- 4 . 64 + 2 . 22 “* 

6.7 

5012 

NGC 6281 

17:04:41 

-37:59:06 

347.731 

1.972 

0.479+0.096 

-5.58+0.26 

-3.43+0.55 

-3.60+0.55 

7.5 

314 

NGC 6404 

17:39:37 

-33:14:48 

355.659 

-1.177 

2.400+0.480 

10.60+1.10 

- 0 . 67 + 2 . 93 ^ 

- 0 . 55 + 3 . 83 “* 

5.6 

501 

NGC 6583 

18:15:49 

-22:08:12 

9.283 

-2.534 

2.040+0.408 

3.00+0.40 

+ 0 . 32 + 3 . 54 ^ 

- 1 . 44 + 3 . 57 '* 

6.0 

1000 

NGC 6633 

18:27:15 

+06:30:30 

36.011 

8.328 

0.376+0.075 

-28.95+0.09 

-0.21+0.31 

-1.60+0.31 

7.7 

426 

IC 4756 

18:39:00 

+05:27:00 

36.381 

5.242 

0.484+0.097 

-25.80+0.60 

-0.14+0.23 

-3.40+0.23 

7.6 

500 

NGC 6791 

19:20:53 

+37:46:18 

69.959 

10.904 

5.853+1.171 

-47.10+0.70 

-0.57+0.13 

-2.45+0.12 

8.2 

4395 

NGC 6819 

19:41:18 

+40:11:12 

73.978 

8.481 

2.360+0.472 

2.34+0.05 

-3.14+1.01 

-3.34+1.01 

7.7 

1493 

NGC 6939 

20:31:30 

+60:39:42 

95.903 

12.305 

1.800+0.360 

-42.00+10.00 

0.86+0.46 

-2.16+0.53 

8.4 

1585 

NGC 7142 

21:45:09 

+65:46:30 

105.347 

9.485 

2.344+0.469 

-50.30+0.30 

1.06+0.51 

-4.43+0.34 

8.9 

6918 

NGC 7160 

21:53:40 

+62:36:12 

104.012 

6.457 

0.789+0.158 

-26.28+3.88 

-1.57+0.51 

-1.60+0.51 

8.2 

19 

NGC 7789 

23:57:24 

+56:42:30 

115.532 

-5.385 

1.795+0.359 

-54.70+1.30 

4.08+0.72 

0 .21+0.68 

8.9 

1412 


Source: 1) Friel et al. (2010), 2) Jacobson et al. (2011), 3) Pancino et al. (2010) and 4) extracted from UCAC4 Catalogue by Zacharias 
et al. (2013). 
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Table 16. The present positions and space velocity components of 69 OCs considered in this study along with their membership 
probabilities. 


Cluster 

X 

y 

z 

UgSR 

VqsR 

WqsR 

probability 

Reference 


(kpc) 

(kpc) 

(kpc) 

(km see” ^ ) 

(km see” ^ ) 

(km see” ^ ) 

percent 





Thin 

disk 





NGC 752 

-8.308±0.062 

0.285±0.057 

-0.180zh0.036 

-4.0zh 2.0 

207.3± 4.4 

-13.0± 3.6 

98.9zh 0.2 

Paper I 

NGC 1817 

-9.909±0.382 

-0.208zh0.042 

-0.447zh0.089 

-42.5zh 3.4 

163.0zhl4.7 

-32.6zhl0.9 

77.8zh22.0 

Paper I 

NGC 2360 

-9.215±0.243 

-I.443zh0.289 

-0.047zh0.009 

-69.2zhl6.1 

252.7zhl3.4 

13.2zhl3.6 

97.5zh 1.8 

Paper I 

NGC 2506 

-10.161it0.432 

-2.636zh0.527 

0.597zh0.119 

-66.8zh 5.3 

172.6zh 2.9 

-9.8zh 7.1 

93.9± 1.5 

Paper I 

NGC 2335 

-9.024ib0.205 

-0.979zh0.196 

-0.029zh0.006 

23.6zh 6.6 

216.Ozh 6.8 

-8.1zhl0.5 

99.1± 0.2 

Paper II 

NGC 2482 

-8.636ih0.127 

-I.182zh0.236 

0.048zh0.010 

-30.9zh 5.4 

201.8zh 2.7 

-12.5zh 5.4 

98.6zh 0.3 

Paper II 

NGC 2251 

-9.236ih0.247 

-0.541zh0.108 

0.002zh0.000 

-8.4zh 1.8 

201.7zh 4.0 

-1.5zh 3.9 

99.Izh 0.1 

Paper II 

NGC 2527 

-8.243ih0.049 

-0.550zh0.110 

0.019zh0.004 

-26.0zh 4.2 

196.7zh 1.8 

8.8zh 1.6 

98.7zh 0.1 

Paper II 

NGC 2539 

-8.790ih0.158 

-1.079zh0.216 

0.263zh0.053 

-13.2zh 1.9 

199.6zh 1.2 

-14.9zh 5.9 

98.6zh 0.4 

Paper II 

NGC 2682 

-8.625ih0.125 

-0.451zh0.090 

0.480zh0.096 

-28.5zh 3.1 

187.5zh 4.2 

-11.7zh 7.6 

98.Ozh 0.6 

Paper II 

NGC 1342 

-8.581ih0.116 

0.271zh0.054 

-0.176zh0.035 

20.5zhl.3 

217.8zh2.6 

1.0zh3.2 

99.2zh0.0 

This paper 

NGC 1662 

-8.404ih0.081 

-0.055zh0.011 

-0.157zh0.031 

24.6zh0.6 

226.0zh0.6 

6.5zhl.2 

99.2zh0.0 

This paper 

NGC 2447 

-8.516ih0.103 

-0.899zh0.180 

0.002zh0.000 

-27.9zh5.6 

221.2zh3.2 

-2.0zh2.5 

99.2zh0.0 

This paper 

NGC 2354 

-8.750ih0.150 

-I.223zh0.245 

-0.171zh0.034 

2.4zh4.3 

192.9zh2.4 

-20.6zh6.5 

97.8zh0.8 

This paper 

NGC 1912 

-9.056ih0.211 

0.142zh0.028 

0.013zh0.003 

6.5zh0.7 

201.9zh4.8 

-7.0zh3.0 

99.0zh0.1 

This paper 

NGC 1245 

-10.477ih0.495 

1.627zh0.325 

-0.466zh0.093 

48.9zh 4.4 

215.2zh 4.9 

-47.1zhl2.9 

78.4zh29.7 

Jacobson et al. (2011) 

NGC 188 

-8.854ih0.171 

1.320zh0.264 

0.648zh0.130 

40.4zh 8.9 

200.4zh 6.6 

-13.2zh 9.3 

97.4zh 0.6 

Friel et al. (2010) 

NGC 2099 

-9.380ih0.276 

0.055zh0.011 

0.075zh0.015 

-0.3zh 0.5 

172.9zhl0.7 

5.6zh 1.9 

97.2zh 1.2 

Pancino et al. (2010) 

NGC 2141 

-ll.813ih0.763 

-1.248zh0.250 

-0.408zh0.082 

O.lzh 3.5 

177.0zhl0.7 

11.5zh 7.3 

97.4zh 1.2 

Jacobson et al. (2009) 

NGC 2158 

-12.006ih0.801 

-0.471zh0.094 

0.125zh0.025 

-8.7zh 2.6 

154.3zhl4.2 

1.8zh 5.4 

92.6zh 5.7 

Smiljanic et al. (2009) 

NGC 2194 

-9.815ih0.363 

-0.566zh0.113 

-0.078zh0.016 

13.3zh 2.9 

191.2zh 8.4 

-14.4zh 7.2 

98.2zh 0.6 

Jacobson et al. (2011) 

Collinder 110 

-9.114zh0.223 

1.850zh0.370 

0.327zh0.065 

-17.6zh 4.2 

191.Izh 4.8 

4.9zh 4.0 

98.6zh 0.2 

Pancino et al. (2010) 

NGC 2355 

-9.731zh0.346 

-0.752zh0.150 

0.395zh0.079 

-21.7zh 4.6 

195.0zhl2.5 

-16.9zhl0.7 

98.2zh 1.1 

Jacobson et al. (2011) 

NGC 2420 

-10.219zh0.444 

-0.729zh0.146 

0.834zh0.167 

-53.8zh 2.2 

162.2zh 8.8 

1.6zh 7.6 

92.9zh 3.0 

Pancino et al. (2010) 

NGC 2477 

-8.364zh0.073 

-I.241zh0.248 

-0.132zh0.026 

-0.4zh 3.7 

220.7zh 1.1 

6.5zh 3.4 

99.2zh 0.0 

Bragaglia et al. (2008) 

NGC 2660 

-8.195zh0.039 

-2.815zh0.563 

-0.149zh0.030 

-117.2zh27.5 

212.Izh 1.9 

7.0zhl0.9 

92.9zh 6.9 

Bragaglia et al. (2008) 

NGC 2324 

-ll.162zh0.632 

-2.096zh0.419 

0.219zh0.044 

-10.6zhl0.6 

175.0zhl6.1 

-40.0zh20.7 

79.8zh40.1 

Bragaglia et al. (2008) 

IC 2602 

-7.946zh0.011 

-0.151zh0.030 

-0.014zh0.003 

-l.Ozh 2.9 

210.4zh 1.1 

6.9zh 0.2 

99.Izh 0.0 

D’Orazi et al. (2009) 

NGC 3532 

-7.834zh0.033 

-0.463zh0.093 

0.012zh0.002 

-14.7zh 5.3 

211.6zh 1.9 

7.8zh 0.7 

99.Izh 0.0 

Smiljanic et al. (2009) 

IC 2714 

-9.181zh0.236 

-0.388zh0.078 

0.444zh0.089 

-2.9zh 4.4 

236.Ozh 1.8 

-17.8zh 6.4 

98.9zh 0.4 

Smiljanic et al. (2009) 

NGC 3960 

-7.238zh0.152 

-I.674zh0.335 

0.199zh0.040 

-52.1zhl0.7 

218.8zh 5.4 

-12.6zh 4.5 

98.6zh 0.3 

Bragaglia et al. (2008) 

NGC 1883 

-ll.711zh0.742 

1.124zh0.225 

0.419zh0.084 

34.4zh2.5 

215.2zh 6.9 

-33.9zh 10.3 

95.5zh4.3 

Jacobson et al. (2009) 

Collinder 261 

-6.851zh0.230 

-1.852zh0.370 

-0.211zh0.042 

-63.6zhl4.0 

216.0zhll.0 

-3.3zh 8.0 

98.6zh 0.4 

Sestito et al. (2008) 

NGC 5822 

-7.269zh0.146 

-0.577zh0.115 

0.058zh0.012 

-42.3zh 5.9 

204.5zh 7.8 

-8.4zh 3.0 

98.7zh 0.2 

Smiljanic et al. (2009) 

NGC 6134 

-7.172zh0.166 

-0.385zh0.077 

-0.003zh0.001 

-20.6zh 2.2 

220.6zh 4.6 

-3.5zh 4.4 

99.2zh 0.0 

Smiljanic et al. (2009) 

NGC 6192 

-6.540zh0.292 

-0.509zh0.102 

0.057zh0.011 

14.2zh 3.6 

261.4zhl0.6 

1.8zh 8.2 

98.8zh 0.5 

Magrini et al. (2010) 

NGC 6281 

-7.532zh0.094 

-0.101zh0.020 

0.016zh0.003 

2.2zh 0.6 

215.4zh 2.5 

8.4zh 1.3 

99.2zh 0.0 

Smiljanic et al. (2009) 

IC 2391 

-7.999zh0.000 

-0.174zh0.035 

-0.021zh0.004 

-2.8zh 2.6 

213.7zh 0.6 

-19.1zh4.9 

99.2zh 0.0 

D’Orazi et al. (2009) 

NGC 6633 

-7.699zh0.060 

0.219zh0.044 

0.054zh0.011 

-11.4zh 0.5 

206.Izh 0.6 

2.4zh 0.6 

99.Izh 0.0 

Smiljanic et al. (2009) 

IC 4756 

-7.612zh0.078 

0.286zh0.057 

0.044zh0.009 

-6.2zh 1.1 

204.3zh 1.2 

1.9zh 0.8 

99.Izh 0.0 

Smiljanic et al. (2009) 

NGC 6404 

-5.607zh0.479 

-0.176zh0.035 

-0.049zh0.010 

20.0zh 3.5 

215.1zh40.9 

10.0zh36.5 

99.1zhl.O 

Magrini et al. (2010) 

NGC 6583 

-5.989zh0.402 

0.329zh0.066 

-0.090zh0.018 

14.2zh 6.1 

215.0zh34.0 

-2.3zh34.3 

99.2zh0.2 

Magrini et al. (2010) 

NGC 6791 

-6.650zh0.270 

3.707zh0.741 

0.760zh0.152 

38.9zh 9.2 

167.8zh 3.0 

-12.4zh 3.3 

91.8zh 2.4 

Boesgaard et al. (2009) 

NGC 6819 

-7.357zh0.129 

2.244zh0.449 

0.348zh0.070 

58.0zhl4.4 

212.3zh 4.6 

19.9zhll.5 

97.7zh 1.6 

Bragaglia et al. (2001) 

NGC 6939 

-8.182zh0.036 

1.749zh0.350 

0.384zh0.077 

24.8zh 4.9 

189.Izh 9.8 

-17.9zh 5.6 

97.6zh 0.9 

Jacobson et al. (2007) 

NGC 7160 

-8.190zh0.038 

0.761zh0.152 

0.089zh0.018 

24.5zh 2.6 

202.Izh 3.8 

3.Izh 2.0 

99.Ozh 0.1 

Monroe et al. (2010) 

NGC 7789 

-8.772zh0.154 

1.612zh0.322 

-0.169zh0.034 

3.2zh 8.3 

160.7zh 4.3 

7.1zh 5.9 

94.6zh 1.3 

Pancino et al. (2010) 

NGC 6253 

-6.633zh0.273 

-0.620zh0.124 

-0.164zh0.033 

-31.4zh7.5 

207.5zhl5.8 

-0.6zhl6.2 

99.Ozh 0.2 

Montalto et al. (2012) 

Ruprecht 4 

-ll.470zh0.694 

-3.140zh0.628 

-0.437zh0.087 

-39.8zh39.3 

209.2zh42.4 

6.1zh53.4 

98.9zhl.2 

Garraro et al. (2007) 

NGC 7142 

-8.614zh0.123 

2.229zh0.446 

0.386zh0.077 

43.3zh6.3 

190.8zh 3.0 

-45.7zhl0.0 

73.4zh25.9 

Jacobson et al. (2007) 

Berkeley 39 

-ll.411zh0.682 

-3.242zh0.648 

0.838zh0.168 

-26.4zh 52.8 

169.2zh 52.0 

-31.8zh 78.4 

86.3zh89.9 

Friel et al. (2010) 

Berkeley 75 

-13.196zhl.039 

-7.259zhl.452 

-I.766zh0.353 

-69.4zhll9.3 

170.7zh87.1 

-23.7zhl42.8 

87.2zhl20.9 

Garraro et al. (2007) 



Thin-Thick disk 




Melotte 66 

-8.750zh0.150 

-4.112zh0.822 

-1.062zh0.212 

-169.0zh45.7 

237.2zh 9.8 

-7.9zhl7.4 

56.4zh61.9 (Thin) 
42.7zh59.6 (Thick) 

Sestito et al. (2008) 

NGC 2425 

-9.899zh0.380 

-2.396zh0.479 

0.176zh0.035 

-68.5zhl0.4 

152.Ozh 8.0 

-37.3zhl6.2 

40.6zh44.3 (Thin) 
58.9zh43.8 (Thick) 

Jacobson et al. (2011) 




Thick 

disk 




NGC 2266 

-ll.314zh0.663 

-0.458zh0.092 

0.608zh0.122 

34.8zh 2.4 

179.0zhl3.8 

-56.0zhl4.8 

69.6zh51.2 

Paper II 

Tombaugh 2 

-ll.639zh0.728 

-4.816zh0.963 

-0.729zh0.146 

-107.6zh50.9 

171.9zh40.8 

-67.8zh73.8 

97.4zh5.9 

Jacobson et al. (2007) 

Berkeley 73 

-15.884zhl.577 

-5.595zhl.ll9 

-1.605zh0.321 

6.5zh81.4 

76.2zhl09.6 

-39.3zhl25.0 

93.8zh28.6 

Garraro et al. (2007) 

Berkeley 32 

-ll.010zh0.602 

-1.603zh0.321 

0.263zh0.053 

-74.9zh 7.1 

150.1zhl3.6 

-47.2zhl9.1 

85.7zh31.2 

Yong et al. (2012) 

NGC 1193 

-ll.742zh0.748 

2.442zh0.488 

-0.963zh0.193 

54.3zh 8.7 

112.3zhl6.5 

-33.4zhl6.4 

93.9zh 8.2 

Friel et al. (2010) 

Berkeley 22 

-13.584zhl.117 

-2.028zh0.406 

-0.844zh0.169 

-39.0zh32.6 

109.0zh89.9 

-69.3zhl07.5 

96.0zhl6.4 

Yong et al. (2012) 

Berkeley 18 

-12.972zh0.994 

1.454zh0.291 

0.455zh0.091 

-21.1zh 7.8 

145.3zhl7.7 

-128.5zh28.4 

85.8zhl9.6 

Yong et al. (2012) 

NGC 2243 

-10.147zh0.429 

-3.655zh0.731 

-I.379zh0.276 

-61.8zh25.2 

175.7zhll.3 

58.1zhl9.7 

84.7zh41.7 

Jacobson et al. (2011) 

Berkeley 20 

-15.348zhl.470 

-3.205zh0.641 

-2.509zh0.502 

25.1zh24.6 

43.5zh43.2 

-38.0zh31.5 

75.2zh43.1 

Sestito et al. (2008) 

Berkeley 21 

-14.149zhl.230 

-0.746zh0.149 

-0.272zh0.054 

33.0zh 4.9 

79.8zh33.0 

-101.9zh26.8 

71.1zh37.0 

Yong et al. (2012) 

Saurer 1 

-18.752zh2.150 

-7.466zhl.493 

1.698zh0.340 

-38.6zhl21.3 

127.4zhl78.6 

27.5zhl88.5 

70.0zh430.0 

Garraro et al. (2004) 




Thick disk-halo 




Berkeley 33 

-12.540zh0.908 

-4.622zh0.924 

-0.524zh0.105 

-142.3zh20.2 

358.4zh42.4 

-88.5zh36.2 

56.4zh55.7 (Thick) 

Garraro et al. (2007) 




Halo 



43.6zh55.7 (Halo) 

Berkeley 29 

-22.000zh2.800 

-4.566zh0.913 

2.077zh0.415 

57.7zh22.6 

-71.2zh71.1 

-143.6zh61.5 

100.Ozh 0.0 

Sestito et al. (2008) 

Berkeley 31 

-15.384zhl.477 

-3.655zh0.731 

0.739zh0.148 

-31.2zh 9.5 

134.6zh22.5 

-207.0zh48.4 

98.5zh 6.6 

Friel et al. (2010) 

Berkeley 25 

-15.706zhl.541 

-8.178zhl.636 

-I.922zh0.384 

0.5zhl26.6 

79.2zhl29.4 

-134.3zh211.7 

66.4zh297.4 

Garraro et al. (2007) 
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